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Abstract 

We pursue applications for symplectic Plancherel growth based on a repulsion phenomenon 
arising in its diffusion limit and on intermediate representation theory underlying its correlation 
structure. Under diffusive scaling, the dynamics converge to interlaced reflecting Brownian 
motions with a wall that achieve Dyson non-colliding dynamics. We exhibit non-degeneracy 
of constraint in this coupled system by deriving a path property that quantifies repulsion 
between particles coinciding in the limit. We then identify consistent series of Plancherel 
measures for Sp{oo) that reflect the odd symplectic groups, despite their non-semisimplicity. 
As an application, we compute the correlation kernel of the growth model and investigate its 
local asymptotics: the incomplete beta kernel emerges in the bulk limit, and new variants of 
the Jacobi and Pearcey kernels arise as edge limits. In particular, we provide further evidence 
for the universality of the 1/4-growth exponent and Pearcey point process in the class of 
anisotropic KPZ with a wall. 
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1 Introduction 

Plancherel growth processes are certain continuous dynamics for Gelfand-Tsetlin patterns that can 
serve as prototypes for studying probabilistic phenomena in systems involving short-range depen¬ 
dence and simple interactions (for us, blocking and pushing). Figure [^depicts an instance of the dy¬ 
namics studied in this paper. In our orthogonal/symplectic patterns, level n G Z>o := {1,2,3...} 






Figure 1: Three possible initial steps of either orthogonal or symplectic Plancherel growth. The 
arrows indicate the direction that the particle will attempt to jump next; see Figurefor the fourth 
step. The tiling determined by their positions suggests a stepped-surface interpretation. 

has rn := [(n -|- l)/2j particles |/f, 1 < i in Z>o := (0,1,2,...}. Symplectic Plancherel growth 
features these particles performing independent simple random walks, but a unit jump of a particle 
yf is suppressed or forced if otherwise it would land outside of the interval (write = oo) 

J [O, y'jPi) , i = Tn with n > 1 odd 
1 Vi-i) , otherwise for 1 < i 

In particular, the “wall-particles” y"^^ at odd levels n > 1 are suppressed by the wall at 0. This 
blocking and pushing of a particle by straddling particles on the previous level (including 0 and cxd) 
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Figure 2: The fourth step illustrates the distinctive wall behaviors: reflection (left) and suppression 
(right). Simulations of the orthogonal and symplectic cases indicate similar average behavior. 


serves to maintain interlacing conditions, which in the chosen coordinates read as 


n+1 
i ’ 


0 < Vi < y, 

vm < y" < 


i = Tn with 77 , > 1 odd 
otherwise for 1 < z < 


( 1 ) 


Orthogonal Plancherel growth studied by Borodin-Kuan [T3] differs only by the reflection of jumps 
into the wall at 0. Although the difference between the two systems may seem negligible at first, 
note that a single wall-particle’s movement can cause the models to diverge by infinitely many parti¬ 
cle positions; see Figure |^for a comparison. Since local limit behavior of such integrable probability 
models can be sensitive to their basic parameters (notably, to initial conditions; see, e.g., [H]), it is 
worth investigating to what degree the distinct wall behavior is felt in the correlation structure and 
in various asymptotic regimes. To explain terminology, when deriving these dynamics via repre¬ 
sentation theory, the wall behavior originates from the initial conditions of orthogonal polynomials 
used to express irreducible/indecomposable characters of even-orthogonal/odd-symplectic groups; 
see Section 5.3 This remark indicates another, a priori signihcant departure from the orthogo¬ 
nal case: the relevance of the representation theory of non-semisimple Lie groups intermediate in 
the classical symplectic series. For readers less familiar with these algebraic items, note that all 
necessary concepts are reviewed and that Section relies only on modern probability theory. 


1.1 Main goals and motivation 

We have two primary goals: 
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1. To quantify repulsion in the Dyson non-colliding dynamics achieved in the diffusion limit in 
order to exhibit non-degeneracy of constraint in this regime. 

2. To compute series of Plancherel measures for Sp{oo) that reflect the intermediate representa¬ 
tion theory of the non-semisimple odd symplectic groups, and subsequently investigate their 
role in the correlation structure and local limit behavior of symplectic Plancherel growth. 

In pursuing the hrst goal, we resolve a basic version of the conjecture posed after Proposition 2 
of [12], conhrm the decompositions (4), (29) in [47] when the interlaced system starts at the origin, 
and thus also lift the restriction in [27] that the limiting system must start in the interior of the 
chamber. To show that the amount of constraint in the pre-limit system does not degenerate in 
this regime, we need to quantify the repulsion between particles coinciding in the limit. The path 
property that assumes this role is interesting in its own right: it delineates a lower envelope at 
t = 0 for a multivariate self-similar process with positive components that do not drift to infinity 
almost surely. For such ]R_|_-valued Markov processes drifting to infinity, see [21] for rather general 
techniques and integral tests. The results of Section 2.4| at the heart of this development are new 
in this literature and are inspired by Burdzy-Kang-Ramanan [20]; we otherwise briskly follow the 
elegant procedures of Gorin-Shkolnikov [27] and Warren m- For the related case of vicious random 
walkers, see the series of results by Katori, Tanemura, et.ah, whose paper [21] covers convergence 
of our even levels. 

For the second goal, the non-semisimplicity of odd symplectic groups is a potentially serious 
obstacle for the identification of suitable “intermediate” Plancherel measures for Sp{oo) - such 
measures can no longer arise as coefficients of irreducible characters in the orthogonal decomposition 
of restricted infinite extreme characters. Circumventing this lack of irreducibility relies on the 
independent works of Proctor [4T] and Shtepin [13], and we suspect there are more connections in 
this vein to be made between probability theory and nonclassical representation theoretic objects. 
Defosseux [23] and Warren-Windridge [IE] have already considered our dynamics for symplectic 
Gelfand-Tsetlin patterns, with the latter paper anticipating our interest, but the results we derive 
here are much in the same spirit as Borodin-Kuan [H], who cover the orthogonal case. Obtaining 
rehned local limit behavior in their manner is well codihed, originating with a law of large numbers 
result for the limit shape of random Young tableaux, derived independently by Vershik-Kerov [IB] 
and Logan-Shepp The former pair were studying asymptotics of the maximum dimension 

of irreducible representations for the symmetric group, and the latter were pursuing a variational 
problem arising in combinatorics posed by Richard Stanley. Baik-Deift-Johansson [3], and shortly 
after Okounkov [3H] , derived more detailed results concerning local behavior at the edge of the limit 
shape, inspiring much research on such behavior, including the present work; for a perfunctory list, 
consider [iliniElElElCniilllliniElIlsliiniElEilEHlETlEHllaHlSQlSe]. 

Important Notational Convention. Quantities with indices that overflow are set to 0 and those 
with indices that underflow are set to oo. 


1.2 Main applications 

1.2.1 Diffusive scaling limit: non-degeneracy of constrained dynamics 


The basic description of symplectic Plancherel growth (considered up to level n > 1) readily leads 
to the candidate diffusive limiting system (W(f))i>o with components satisfying 


Wf is a Brownian motion reflected in the interval 


[0,oo), 


if fc = 1 
if fc > 2 


( 2 ) 
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for 1 < /c < n, 1 < z < Tfc, where the are independent and started at 0; note our notational 
convention puts additional static particles at 0 (for index overflow) and at oo (for index underflow). 
Section 2d^ makes the notion of reflecting Brownian movements in a time-dependent domain rigor¬ 
ous, but we say the coupled dynamics of (W(f))t>o are non-degenerate if its components uniquely 
solve the system of equations 


B! + otherwise for 1 < z < ’ 


(3) 


where Ly is twice the semimartingale local time at 0 of a semimartingale Y and where = 0 
(see Remark 4.2 of m for an explanation of the appropriate local time scalings). Intuitively, these 
equations indicate the amount of constraint in the system (kV(t))o<t<T is hnite for every T > 0. But 
demonstrating non-degeneracy of (kV(t))i>o requires that we resolve a subtle yet concrete difficulty. 



Figure 3: The hrst path in black is a standard brownian motion reflected at 0. Reflecting the blue 
path to stay between 0 and the black path requires an “inhnite amount” of constraining to achieve. 



Figure 4: Diffusion limit of hrst three particles nearest the wall (visualized in the single state space 
M+). The level 1 particle in red creates a repulsion between the wall at 0 and the level 2 particle in 
black; although this extra room may still not be enough, we show the rehected level 3 wall-particle 
in blue has non-degenerate constrained dynamics. 

To see the issue in our case, write the dynamics of a typical rehected particle in the basic 
decomposition W = B Y, where the constraining process Y “pushes” the Brownian particle 
B to remain between some independent particles £, r on the previous level. If ever the interval 
[£{t),r{t)] collapses, Y might need to “push an inhnite amount” in order to keep B at that point. 
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For example, Y can be finite variation on every period between collapse times of [i{t),r{t)], bnt 
accnmulate infinite variation over infinitely many snch periods in finite time; see Figure]^ for snch 
a scenario. 

For a more snrprising and relevant example, the interval [£{t),r{t)] can “close too qnickly” at 
t = 0, say, for some c > 1, almost snrely 


lim snp 

t4,o 


r{t) - i{t) 


0 ; 


(4) 


in this sitnation, even if the particles i{t),r{t) remain apart for t > 0, Y will not be of bonnded 
variation. If Y is of infinite variation, then the symbols indicating local times in ([^ do not make 
sense (even thongh their increments may). Proposition 4.13 and the proof of Proposition 4.12 of 
[20] provide more details on these two degenerate examples. 

Fortnnately, by the construction of our system (|^, any particles i,r that can collide must 
straddle an independent Brownian particle on a further previous level that could create a sufficiently 
strong repulsion that Y only needs to “push a finite amount” to keep B between them; see Figure 
Our strategy, then, consists of quantifying this degree of repulsion in order to rule out the 
degenerate possibility (|^ and exhibit non-degeneracy of the constrained dynamics (W(t))t>o of 
(|^. More broadly. Section will establish the following result. 

Theorem 1.1. Fix n>l. The interlaced system (W(t))i>o of reflecting Brownian motions with a 
wall defined by ([^ has levels W^, 1 < k < n, that achieve the Dyson non-colliding dynamics 

{ dBi + + Y^ j^i ( fc flyk ) dt, if k even 

‘P' . , (5) 

dBi + \lprk)dL\flt)+ Y. jp [Ykpp + dt, if k odd 


for 1 < i < r^, where the Bi are independent standard Brownian motions. Moreover, each particle 
2<k<n, l<i< rk+i is reflected in an interval satisfying the almost 


sure path property 


lim inf 

qo 




= oo 


for any c > 1/2; (W(t))t>o is then verifiably the unigue, strong solution to ([^ when started at the 
origin. Hence, under diffusive space-time scaling, symplectic Plancherel growth considered up to 
level n converges in law to the non-degenerate system (|^. 


The path property above is readily seen to hold in Warren’s setting m of interlaced reflecting 
Brownian motions (without a wall), which achieve classical Dyson dynamics 


dXl = dBl 


K 

+ 2 


E 


dt 


Aj-Ar 


l<i < N, 


( 6 ) 


where the 5* are independent. In particular, setting k = 2 and arguing as in Section 2.4 confirms the 


decompositions (4), (29) in |17] when the system starts at the origin and thus also lifts the restriction 
in 123 that the limiting system must start in the interior of the chamber {xi > ... > xn}- Note the 
path property for ([^ is compelling only for > 3; in the case N = 2, the single gap is proportional 
to a 3-dimensional Bessel process and the result is classical. 
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1.2.2 Determinantal correlation kernel and local limit behavior 


As we saw in the last section, the diffusion limit primarily cares inductively about interactions of 
three successive levels at a time. To capture the effects of correlations between particles farther 
apart in the system, we now view symplectic Plancherel growth as determining a point configuration 
A’( 7 ); 7 > 0 , in Z>o x Z>o, he., as a random element of A fundamental theorem about 

this point process is that it is determinantal with an explicitly computable kernel when started from 
the “leftmost” initial condition (see the hrst image of Figure [^. To state this result and others that 
follow, let Jk,i/ 2 {x), Jk-i/ 2 {x), k > 0 denote the Chebyshev polynomials of the second and third 
kind, respectively, (see Section 3.1) and write 


Oik ■ = 


1 / 2 , if fc even 
— 1 / 2 , if k odd 


Consider the recoordinatization X ( 7 ) dehned to have components := X^{'y)+rn — {n — k + 1), 

n > 1, 1 < k < Tn (see Section [4T for explanation and visualization). The next theorem is a 
consequence of the similar statement Theorem 4T concerning the general symplectic Plancherel 
point process. 

Theorem 1.2. The kth correlation functions {'f>l}k>i of are determinantal: for zi,..., Zk & 

Z>o X Z>oj ^ 

^(zi, ...,Zk):= P(T(7) D {zi ,..., Zk}) = det[K^{zi, 

where the (nonsymmetric) kernel K'^ is given explicitly, for {s,n), {t,m) G Z>o x Z>o, by 


K'f{{s,n), {t,m)) = 


2 «n + l /2 ^ 

71 27ri 


'-1 


^'yx 


(1 -2,)r„+a„(i^^)l/2 


yy^Js,an{x)Jt,am{'o)' \r 7 ^ 

(1 — uy"'[x — u) 


dudx 


2an+ll‘2 

+ l(n>m)-7- / (1 + xf^'^dx. 


(7) 


TT 




Here, the u contour is a positively oriented (i.e., counterclockwise) simple loop around [—1,1]. 

Although the correlation kernel ([^ of Theorem 1.2 may appear unwieldy at hrst, its double 
integral expression allows for asymptotic expansions leading to rehned limiting results. Dehne the 
region Qf nontrivial hydrodynamic limiting behavior as the collection of locations and times 

with a positive chance of either seeing a particle or not: 


p^iiquid I u,r]) eR+xR+xR+:0< lim [Np])) < 1 


N^OO 


( 8 ) 


Under this hydrodynamic scaling, the method of steepest descent applied directly to the kernel ([^ 
shows that is exactly the region of (r, z/, rj) G Mi where the function 


S{z) = := T 


z + z 


-1 


U log Z + 7] log 


z + z 


-1 


(9) 


has a unique critical point zo = zo(t, z /, 7 ) in El \ D. Similarly dehne to be the region where 

we are sure to hnd particles in this limit, and ^vhere there can be none; more precisely, these 

are the regions where the limit in ([^ is 1 and 0, respectively. See Figure]^ for a depiction of these 
regions. To identify the nontrivial limit in the liquid region, dehne the incomplete beta kernel by 

B(t, no :=^ 
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Figure 5: The hydrodynamic limit capturing average behavior indicates a densely-packed region of 
inactivity (“frozen”) and a central region of activity (“liquid”), with the remaining area unreached 
(“empty”). Notably, the orthogonal and symplectic cases feature the same curves delineating these 
regions, which are given analytically by Proposition 


6.2 


compare with Figure 4 of 


where the contour crosses (0,1) if A; > 0 and (—oo, 0) if k < 0. This kernel serves as a generalization 
of the sine kernel to describe bulk limits of the Schur process in |3D] by Okounkov-Reshetikhin. 

Theorem 1.3. (Bulk Limit: Incomplete Beta Kernel) Let (r, i/, r/) G and let zq = zo{T,u,ri) 

be the unique critical point of m ]HI\D. Assume 7 > 0, {si,ni) G Z>o x Z>o, 1 < ^ < A:, depend 
on N in such a way that 7 ~ Nr > 0, ~ Nu > 0, and r„. ~ Nr] > 0. Assume the differences 

Si — Sj, Ui — Uj are of constant order. Then as N ^ 00 , 

p^((si,ni),..., (sfc,nfc)) -)■ det[B{ni - np, (si - Sj) + (n* - Uj) - {r^ - rnj)\zo)]lj=i. 

At the edges of the liquid region -we hnd variants of the Jacobi and Pearcey kernels of 

Borodin-Kuan The hrst of these types arises in the large time limit at a hnite distance 

from the wall, away from the corner of the phase transition; it only seems to have appeared so far 
in the similar contexts of na and 135 , though with different Jacobi polynomials. 

Theorem 1.4. (Edge Limit: Discrete Jacobi Kernel) Assume'y > 0, (sj,nj) G Z>oxZ>o; 1 < * < A:, 
depend on N in such a way that 7 ~ Nr > 0, r„. ~ Nrj > 0 but the Sj are fixed and finite. Assume 
only the differences n* — Uj are of constant order. Then as N ^ 00 , 


^((si,ni),...,(sfc,nfc)) ^ 


det[L{{si,ni), {sj,nj)\l - 

1 , 


*/ 1 - v/r > -1 
*/ 1 - v/^ < -1 


where the discrete Jacobi kernel L is given explicitly, for {s,n), {t,m) G Z>o x Z>o, by 

L{{s,n),{t,m)\e) : = - / [l]n>m) + + x^^'^dx. 

^ J-1 

( 10 ) 

Our last result shows the corner of the phase transition (i.e., near where the regions and 

'phquid ] 2 ieet the wall) exhibits a 1/4-exponent in its fluctuation scale, providing additional rigorous 
evidence of the universality of this growth exponent in the class of anisotropic KPZ with a wall. The 
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evolving stepped-surface interpretation indicated by Fignre [fallows one to define a height fnnction 
as “the nnmber of particles at that point and to the right”. The fiuctnations abont its mean should 
be described (at least formally) by the anisotropic KPZ equation, a stochastic heat equation of the 
form 

dth = u^dlh + Uydyh + ^Xxidxh)‘^ + ^Xyidyhf + 

where is a space-time white noise and A^;, Xy have different signs (should they have the same 
sign, this equation reduces to the classical KPZ equation). We do not explore this equation in 
this work, we only note its relevance for future research; see, e.g., Borodin-Ferrari [10] for more 
discussion. Such rigorous results for systems with a wall still seem to be rare; the author only knows 
of [HI [m EH ES] • The kernel arising in our setting at this critical corner point appears to be an 
instance of a general family of Pearcey-type posited in Remark 1.5 of [55| . 

Theorem 1.5. (Edge Limit: Pearcey Kernel) Let denote the kth correlation function of the 
determinantal process determined by the complementary process Assume 7 > 0 , {si,ni) G 

Z>o X Z>o, 1 < i < fc, depend on N in such a way that 7 ~ N/2, s i i > 0, and rn^ — N 

y/Nrii > 0. Then as N ^ 00 , 

^k/4 . ^ (sfc,nfc)) ^ det[/C((z/i, 7 i), (z^j, 


where the Pearcey Kernel K. is given explicitly by 


(^^ 2 , 72 )) : = l(r,i>772) 


V^(7i -V 2 ) 

V 2 1 


exp 


(z/i + z/ 2 )^ 
Vl - V2 


— exp 


(Z/I - t/ 2 )' 

Vl - V2 


TT 2ni 


loo POO • r /o~7i dx'du' 

/ e 8 2 sinyuiy 2 x'\sm.[v 2 v 2u \- 


' —200 J 0 


\/v({u' — x') 


Acknowledgements. The author would like to thank Rene Carmona, Allen Knutson, Daniel 
Lacker, Ramon Van Handel, and Mykhaylo Shkolnikov for many helpful discussions, as well as Jon 
Warren for an informative private correspondence. I would especially like to thank Jeffrey Kuan, 
for early encouragement and consistently valuable guidance, as well as for supporting a simulation 
of the symplectic growth process at http://www.math.harvard.edu/~jkuan/Cerenzia.html, 
Part of this research was performed while the author was visiting the Institute for Pure and 
Applied Mathematics (IPAM), which is supported by the National Science Foundation; the author 
was also partially supported by NSF; DMS-0806591. 


2 Diffusive scaling limit and repulsion of Dyson gaps 

2.1 Time-dependent Skorokhod problem and Skorokhod reflection map 

The machinery of the extended Skorokhod map provides a remarkably natural and careful con¬ 
struction of our particle system and its diffusion limit. A solution to the Skorokhod problem (SP) of 
reflecting a particle in a (potentially time-dependent) interval [£, r], £ < r, is a pair ( 0 , rj) where v 
“pushes” the particle 0 to remain in [i, r] and the resulting reflected movements 0 = 0 J- G [i,r] 
should be regular; classically, v should be of bounded variation and ( 0 , v) should satisfy the familiar 
conditions 

poo poo 

j ^{l!(s)<<f>(s))dv£{s) / l((/i(s)<r(s))d7r (s) 0, (H) 

Jo Jo 
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where are nondecreasing functions affording the (minimal) decomposition rj = rji — rjr. But as 
explained in Section 1.2.1, the conditions 0 are too restrictive for our situation. To state a weaker 
notion of solution, let Z1(]R+, S) be the space of right-continuous functions with left limits on IR+ that 
take values in a separable complete metric space S. Write D := D{R+, M), D~ := D{R+, [— 00 , cso)), 
and := D(K^, (— 00 , 00 ]). 


Definition 2.1. (Dehnition 2.2, [20]) Given {i, r^ip) E D xD^xD with i < r, a pair (</>, 77 ) G D x D 
solves the extended Skorokhod problem (ESP) on [£, r] for -0 if 


1 . 0 (f) = ' 0 (f) -h 77 (f) e [i{t),r{t)] for f > 0 

2. For 0 < s < t < 00 , 

77 (f) — 77 ( 5 ) > 0 , if 0 < r on (s,f]; and 77 ( 5 ) — 77 ( 5 —) > 0 , if 0 (s) < r(s) 
77 (f) — 77 ( 5 ) < 0 , if £ < 0 on (s,f]; and 77 ( 5 ) — 77 ( 5 —) < 0 , if i{s) < (j){s) 


Given {£, r, 0) G D x x D with £ < r, an explicit candidate for the reflected process 0 is 
provided by the extended Skorokhod map r(f', r|0) := 0 — Ee^ri'ip), where 


S£,r(0)(f) := max I (0(0) — r(O))’'' A inf (0 (m) — £{u)), sup [(0(s) — r(s)) A inf (0 (m) — £{u))\ 

\ «e[o,t] s&[o,t] u&[s,t] 

Note some simplihcations: if r = cx), then S£^oo(0)(0 = ~ supgg[o, 7 ][^(' 5 ) ~0(s)]'''; ^ > 0; also, at the 
initial time 0 , 

S,,00)(O) = [0(0) - r(0)]+ A [0(0) - 00)]-. (12) 

In words then, r(£, r|0) starts by putting the particle 0 at the nearest point in [£,r] and maintains 
its increments except for putting those that land outside of [£, r] again at the nearest point inside. 
Theorem 2.6 and Proposition 2.8 of |20| show P is a jointly continuous function on D~ x x D 
(each topologized by uniform convergence on compacts) and provides the unique solution 

( 0 , 77 ) := (r(0r|0),r(0r|0) -0) 


to the associated ESP on [£,0 for 0. The condition in0g[t^^t^0r(f) — £{t)) > 0, 0 < fi < f 2 , is 
sufficient for ( 0 , 77 ) to satisfy the more regular conditions (11) on [fi,f 2 ]; cf. Gorollary 2.4 of [20] . 


2.2 First construction of symplectic Plancherel growth and continuum 
universal limit 

Fix 77 > 1. Gonsider the collection X{t) := {Xj^(f)|l < k <n, 1 < i < r^}, f > 0, of continuous time 
simple random walks X!^ started at 0, with jumps dictated by independent unit rate poisson pro¬ 
cesses. For a parameter iV > 0, consider the diffusively scaled versions X(f;iV) := XtN/\fN- Gon- 
struct a process (T’(f, N))t>Q with components X^{t, N) driven by X by setting := P(0, cxd|X0 
and inductively for 2 < k < n, 

fr(O,A;0-0 - 1/0]V|X^), if i = rk with k odd 

* ■ + l/\/iV, A’0-0 - l/y/N\X^), otherwise for 1 < f ^ ^ 

where our convention = 00 applies. This procedure determines deterministic maps {X){t) = 
X{t,N). 
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Definition 2.2. Symplectic Plancherel growth considered up to level n > 1 is characterized by the 


interlaced dynamics of the process X{pf) := ^’( 7 ,1), 7 > 0, defined by (13). 


Note 2.1. By the simplification (|I^ at time t = 0, df(0) is exactly the “leftmost” initial condition 


up to level n. However, Definition 2.2 refers to the interlaced dynamics and does not depend on 


this initial condition. In fact, the paper [27] illustrates how the results of this section hold for much 
more general dynamics and initial conditions, and a forthcoming technical report will show the 


dynamics of m EaiEiiisi m] all converge to the same system under diffusive scaling. It is only 
for continuity of focus and simplicity of presentation that we maintain our concrete setting. 

The appropriate continuum state space for the diffusion limit is given by the collection paths 


of arrays {y^ G M+|l < k < n, 1 < i < r^} of nonnegative real numbers satisfying 

< vi < !/r‘. !<*:<»■«. 


Consider the collection B{t) := {Bf{t)\l < k < n, 1 < i < r^}, t > 0, of independent Brownian 
motions started at 0. Construct now a process (W(f))t>o in paths with components Wf driven by 
{B{t))t>o by setting >V| := r( 0 ,oo|i?;[) and inductively 



(14) 


for 2<k<n, l<i<rk (recall our notational convention). This procedure determines another 


deterministic map := W. Notice as —>■ 00 . 

Proposition 2.3. Fix n > 1 and write dn '■= rn{rn + l)/2 + r„_i(r„_i + l)/2. Then considered 
up to level n, the diffusively scaled symplectic Plancherel growth process {X(t-,N))t>Q converges in 


law on D([0, 00 ), M*^") to the interlaced reflecting system (W(f))t>o of (14) with a wall started at 



Proof. For n = 1, the laws of Xl converge on Z1([0, 00 ), M) to the law of a standard Brownian motion 
(see, e.g., the much stronger Theorem 16.14 of [22], where convergence in the Skorokhod topology 
strengthens to convergence in the supremem norm topology since the limit is continuous). In turn. 


the law of the images X{flN) = r( 0 , cxd|X|(-, iV))(t) converge to a standard reflected Brownian 
motion f ( 0 , oo|Hj)(t) = ^^^(t) + sup,g[o,t] (s)]^- 


For n > 2, assume that the statement of the theorem holds for all levels fc, 1 < fc < n — 1, so 
that the restriction of X{t]N) to these levels converges in law on D([ 0 , cxd), As iV —)■ 00 , 

the rescaled converge in law on D([ 0 , 00 ), M), in the supremum norm topology, to Bi, 

1 < i < Vn, where the Bi are independent standard Brownian motions started at 0. Independence 
ensures the component-wise topology of uniform convergence on compacts is respected. Since 
the Skorokhod map is jointly continuous in this topology (cf. again Theorem 2.6 of [20]), the 
continuous mapping theorem in the form of Theorem 4.27 of [29] implies that as —>■ 00 the laws 
of {Xfl{t] N))t>o, 1 < i < r„, converge on D([ 0 , cx)),M) to (VV"(t))t>o, which completes the proof. 


□ 


2.3 Dyson non-colliding dynamics of limiting system levels 


To determine the dynamics of the limiting system (W(t))i>o constructed in (14), we follow Warren’s 


approach ra of working with two levels at a time. Keep n > 1 fixed and fix fc, 1 < fc < n. Define 


:= {u eMff : Ui > ■ ■ ■ > >k 0 }, 
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(15) 






For u G and v G let m -< n denote interlacement: 


> Ml > ■ • ■ > 't’rfe > Mrfc > Vr^^+1, 

and additionally define 

y^k+i,k ^ y^fc+1 :u^v}, (16) 

For fixed {v,u) G PF^+i’^, consider a filtered probability space (p‘,J^,{J^t}t>o,Q^yu)^ snpporting 

a pair {Vt,Ut) of x M(j?-valned interlaced processes, i.e., f/t -< V) for t > 0, with the follow¬ 
ing dynamics nnder Take /3i,...,to be independent J^t-Brownian motions started at 

Ml,..., Let Uiit) := \l3iit A r)| to be stopped at the collision time 

r := inf{t > 0 : |A-i(f)| = |A(^)I, for some i, 1 < i < r^+i}. (17) 

Then take 7i,... ,7rj.+i be independent J^-Brownian motions started at mi, ... and let the 

dynamics of V be reflected on U throngh the extended Skorokhod map by 

Vi{t) := r{Ui, Ui-il'jii- A r))(f), 1 < i < Vk+i- 


(note onr notational convention applies thronghont the above definitions). Since U starts at strictly 
separated components u G hF^, V admits, for t < r, the decomposition 


^ - ^(’c/i_i-y)(^)’ * = ^k+i with fc + 1 odd 

+ ^v,_uy{t) - ^^._^_yy{t), otherwise for 1 < i < r^+i ’ 


(18) 


(see Remark 4.2 of [20]; cf. also [12] )• These expressions allow us to derive the transition probabilities 
for {V,U) on explicitly; write := ^and 4>t(a:) := J^^(ptiu)du. 

Proposition 2.4. Fort > 0 and {v,u), {v',u') G transition density 

Q{v,u){^t e dv',Ut G du,t < t) = q^{{v,u), {v',u))dv'du' 

for the process (F, U) killed at time r is given by 


Qhi'v,u), (m',m')) := det 


A^{v,v') B^{v,u') 

cf{u,v') Df{u,u') 


where 


' {A’l{v,v'))ij = - Vi) F (-l)^0t(M' + Vi), I <i,j < Tk+l 

{B^{v, u'))ij = $i(m) - Vi) -1- (-l)*^<Fi(M) Vi) - lo+i<i) I <i < Vk+i, 1 < j < 
{C^{u,F))ij = - Ui) + (-l)^"Vfn^i + “*) ^ <i <rk: I < j < Tk+l 

XD^{u,u'))ij = - Ui) + {-l)^-^(j)t{Uj + Ui), l<ij <rk 


The presence of the indicator “l(j+i<i)” in the definition of B^ can in part be explained by the 


calculation required for (23). 


Proof. The proof follows the same line of reasoning as Proposition 2 of m- Fix w' = {v',u') G 
^rk+l,k define G{t,{v,u)) = iv,u)) := q^{{v,u),{v' ,u')). We need to check that G 

satisfies the heat equation on (0, oo) x x and the appropriate boundary conditions. We 
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focus on the latter; the former follows since all functions of {v, u) involved in the expression above 
solve the heat equation. 

In the reflecting wall case k odd, we compute \ur^^=o{Df{u, u'))r^i = (u, n'))rfcZ = 

0, for any 1 < / < so that du^Jur^=oG(t, {v,u)) = 0. For the absorbing wall case k even, 
we have G(t,{v,u)) = 0 when Ur^. = 0. For either parity of k, ii Vi = Ui, 1 < i < r^, we have 
d^^\y^=uiA{v,v')ij = -Gt{u,v')ij and dAvi=uiBt{v,u')ij = -Dt{u,u')ij, so that iv,u)) = 

0. The case Uj+i = Ut is the same. Now if Mi_i = Ui, 1 < i < r^, then the u, and Uj+i rows of 
[G^{u,v'),D^{u,u')] are equal, giving G{t, {v,u)) = 0. Finally, one readily checks the correct initial 
condition lim4o G'(^, {v,u)) = 

Now £x / : —>■ M bounded and continuous. Note that the function 

F{t,w) := [ q^{w,w')f{w')dw'. 

Jy/k+l,k 

solves the heat equation on (0, cxd) x and satishes the same boundary conditions as above. 

Hence, for T, e > 0, Ito’s Formula shows that the process (F(T—f+e, (Vt, Ut)), Ft)t&\o^T] is a bounded 
local martingale, and thus a true martingale. This property, the bounded convergence theorem, 
and the regularity of (cf. Lemma 1 of |17]) together give 


/ qxA^ w') f {w')dw' = limF(T + e, (n, u)) 

JlYk+i,k 40 

= limEgfe [F(e, (Ft, tfr))l(T<T)] = Egfc [/(Fp, t/r)l(T<T)] 

completing the proof. □ 

Let be the law of the stopp ed p rocess Ut started at n G IF^ with the components Ui{t) = 
\/3i{t A r)| from above. Proposition 2.4 implies that for u' G 

Pu{Ut G du',t < t) = p^{u,u')du := det[0t(M'- - Ui) + . 

Now define 


hk{u) := 


ni<i<j<rfe(“?-“i) ■ if k even ^ J 1/2, if k even 


n 


l<i<j<rfe 




if k odd 




— 1/2, if k odd ’ 


(19) 


Note the functions A; > 1, are positive harmonic in IF^ with respect to the generator of Ut, 
and u G IF^ has strictly separated components by definition. We may then define the Doob hk- 
transform P^’+ of (cf. [28]) by 


pk,+ hk{U{t)) _ 


hk{u) 


with hfc-transformed density 


Bu~^{Ut G du',t < t) = Pt’~^{u,u')du' : = 


hkW) 

hk{u) 


Pt{u, u')du'. 


( 20 ) 


Under this measure, r is inhnite and expression (20) implies, by Lemma 3 of [32| and Ito’s formula 
for convex functions (Theorem 7.1.(v), [301), that Ut satishes the non-colliding dynamics ([^, cor¬ 
responding to the case k. We further use Lemma 4 of [32| (cf. also |3]) to arrive at the law Pq’^ of 
Ut issued from the origin: 


Po''’+(Pi G du') = At{u')du' : = 


l-rkiTk+otk) 


Gk 


exp 


\u 


/|2 


2t 


■ hkiu'Ydu' 


( 21 ) 
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Now for w = [v, u) G the marginal u G has strictly separated components, so we 

may define a measure to be the Doob hfc-transform of Q^\ 


Q 


k,- 


hkiUit)) , 
hk{u) 


on 


with hfc-transformed density 


e dw') = q^'^{w,w')dw' := ^^^^\ q^{w,w')dw' 

hk[u) 


for w' = {v',u') G f > 0. Note for v G 


' {u^W^:u^v^ 


hk{u)du = 


hk+i{v) 

k\\ 


so that A^(u,m) := is a transition kernel from down to W^] in particular, A^(u, •) is 

a probability density on {u G : u -< u}. This density together with the entrance laws 
furnish an entrance law k'tdv, u)) := /i^'''^(u)A^(u, u) of a measure Q^q~^q) for the process (V, U) issued 
from the origin: 


^ = ^tiw)dw, z/f (w) = / u^{w')qt’^{w',w)dw' 


'(0,0)^ 


llYk+l,k 


The proof that iy^{{v,u)) = is an entrance law is a consequence of an intertwining 

relation that holds between the families {p\’^}k>i and {qt'""]k>i, t > 0: 


I X^{v,u)q^'~'~{{v,u),{v',u))du'= {v',u)eW 

J {u^W^ :u^v} 

( 22 ) 

This intertwining readily follows by directly “integrating u out” in the expression for q^ in Propo¬ 
sition [231 

[ (lti{v,u),{v',u))du = p^t+\v,v'). (23) 

J {u<^W^ \u^v} 

We now use these properties to determine the projected dynamics of V. 

Proposition 2.5. Under marginal {Vt)t>o has distribution , so the process V) 

satisfies the system ([^, corresponding to the case (A; -f 1). 


/ /\\ J. f _ r, ' 


rk-\-l^k 


Proof. Take 0 < fi < ■ ■ • < tm- Then iteratively using the intertwining relation (22), compute 


^?o'o'|(^ii ^ ^1’ ^*2 ^ ^2, • • • , Vt^ G Am) 


h0,0) 


= / dvi dui--- dVm dUm- Jyti{'^l,Ui)Y\_qti\_^{{p-l,Ui-i),{Vi,Ui)) 

J Ai J {u\^v\} JA 771 


i=2 


= / dvi 
Jai 


P m « 

I dVm ' Pti (^l) J_ J_ Pti—ti-i (^i—1) "^i) ■ / {'VmiUm)dUm- 

-^m 2 — 2 


Noting A^(um, •) is a probability density on {um G : Um -< Vm] completes the proof. 


□ 
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2.4 A path property of Dyson gaps and non-degeneracy of coupled 
dynamics 


Proposition |2.5| does not explain how the process Ut is involved in the non-colliding dynamics of Vt- 
In this section, we conhrm that the non-degenerate constraining dynamics (18) coupling Vt with 
Ut survive sending the starting point to the boundary of the chamber; thus, under the given setup 

the coupled system (Vt, Ut) is a genuine semimartingale. 

Besides being realizable as a Doob transform, Ut is also the eigenvalue process of a certain 
matrix-valued process iff; see Section IV of [T] or Section III.B of [32] • In particular, the process 

inherits the Brownian scaling property Ui{t) = \/tUi{l) and we may cite the Hoffman-Weilandt 
inequality (Lemma 2.1.19, [2]): 


' K 

y ([/,(() - Ui(s) f < 


(24) 


i=l 


Proposition 2.6. Fix k > 1, 1 < i < and h < 1. Write 51 := Ui_i{t) — Ui{t). Then for 


z < e 


P ( inf 5i <z] < 


(log 2 )^ 


Proof. For x in the chamber of (15), write 






+ 


rr* , rr* , I rr* , 

tL/2 J i J 


+ 


1(A: even) 


U 




^2{x) := -d^^Ui^ifx) > 0. (25) 


Let \{f) = (Ai(t) > ... > Kk(f)) denote the unique strong solution of 

d\i{t) = dl3i{t) -f Ui^i{X{t))dt, I <i <rk, 

started at the origin, and notice A agrees with the process U under Pf’"'", except that Ur^ = |Aj,j.| 
for k odd. Note the identities (see pg. 252 of 0) 


rk / \ rk 

'^Ui^l{x)Xi = 21 M +rkl)k even), '^[W2{x) - Ui^i{x)‘^] = 0. 
i=l ^ ' i=\ 

Consider the Lyapunov function 

2 2 r 

fix) ^ 0)2 _ _ ^ [log(a:i - x^) F log(a:i -f Xj)] -V l(fc ei>en) ^ log 

i=l *=1 

and note the inequality (see pg. 251 of Bl) 

< g^i(/0+8) fQj, 1 < i < j < Tfc+I, 


Xi 


(26) 


Xi ± Xj 

where our notational convention applies, i.e., Xr^+i = 0 if fc even. We can also readily calculate 


(27) 


dx fix) = — 

* Tfc 


i 

2xi - —Ui,i{x) 


> dlfi^) = - 


Tk 


i 

2 H-Mi,2 (a:) 

Tk 
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Putting these facts together with Ito’s formula gives for t E [l,h 


rk 


rk 


rt r 




k Jl 


2=1 


2Xi{s) -Mi,i(A(s)) 

rk 


d/3i{s). 
(28) 

where Nk := (2rfc - 1) + 2 ■ l(fc even)- In particular, Zt := /(A(f)) + ^ ^i, 2 (^(s))ds is an 

K ^ 

J^rsubmartingale for f G [1, h~^]. Hence, using (27), we have for 

E[/(Ai) + 2iVfc-(h-i-l) + M(h-i)]' 


P ( inf (5i < z 1 < P 


sup Zt> \log(- ] - 


l<t<h- 


< 


3 log ^ + 


where M{t) is the martingale term of (28) and where we have used Doob’s submartingale tail 
inequality with parameter p = 2 conditional on J^i. To see the expectation is hnite, hrst note the 
term E/(Ai)^ is hnite by applying the inequality log(l/a;) < a(l/a;)^/“ for a, a; > 0: for example, we 
can estimate 


E(loghl)^ = E[(log(l/hl))^l(,j<i)] + E[(log(hl))^l(,j>i)] < 4 E 


+ 


< Cfc < CX), 


where the upper bound C(, is independent of t since [1, h is compact and bounded away from 0. 
Moreover, by Ito’s isometry and independence of the J^j-Brownian motions (3i, we have 


rk ph- 


EM{h 


- 1\2 


E 


E 


k j=i ai 


2Ai(s)-ni,i(A(s)) 

rk 


ds < oo, 


which is hnite by the form (21) of the density and by dehnition (25) of the function Ui^i{x). This 
completes the proof. □ 

Note 2.2. To establish this proposition in the previous version of this paper, we hrst used the 
integrability of E[(5*(5J]“'^, q < 3/2, along with basic chaining (cf. [15]) to prove that (log5^)^ is 
Holder continuous on [1, h~^] of parameter /d G (0,1), or when applied to the classical system 0. 
of parameter (3 G (0,1(1 — ;;;qq))- This earlier approach thus fails for k = 1, but our new method 
captures a special probabilistic structure of the Dyson gaps with the Lyapunov function (26) and 
(with some extra ehort) is sharp enough to be adapted to the critical GOE case. 

Proof of Theorem We proceed by hrst demonstrating the path property of the consecutive 
Dyson gaps 5) := f/j_i(f) — Ui{t) and then exploiting it to deduce non-degeneracy of the constrained 
dynamics (18) for V. Fix 1 < i for k >2, and hx G > 0. Then 

P 


( lim inf — < G I = P 
\ 40 P 


inf ^ < 


G, i.o. n ) = P (limsup | inf — < G 

/ V [hr<t<hr-i p 


It is sufficient to show the probabilities of these distinguished events are summable, so we have 


P 


inf ^ < G 1 = P 
h”<t<h"-^ P 


inf ^ < G(h”)"-B 2 | < p 


inf 5i < 


where we have used Brownian scaling for each hxed t in the inhmum and implicitly the fact 61 > 0. 
By Proposition 2.6 and the fact < 1, we may invoke Borel-Cantelli to conclude the path 

property. 
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Now, to reduce notational burden, write W := Vi, B := 7 *, i{t) := Ui{t), r{t) := and 

consider the decomposition W = B + Y. The non-colliding dynamics (|^ satisfy inftg[tj^t 2 ](r(t) — 
£(t)) > 0 for 0 < < t 2 , so the variation V[t^,t 2 ]F of Y on any such interval [^ 1 ,^ 2 ] is hnite 

by Corollary 2.4 of jJO]. It therefore suffices to check the same is true for some interval [ 0 , 7 ), 
rj G (0,1), at 0. We have just shown that for any c G (1/2,1), almost surely 


40 V 


= 00 . 


(29) 


As a result, f{t) := r{t) — £{t) satisfies f{t) > g{t) := V for small f > 0. Following the parabolic box 
approximation of Theorem 4.8 of [20], let {sfc}fcgz<o enumerate, decreasing as k I — 00 , the family of 
points {2~P + q ■ indexed by g = 0,..., [2~p/ g^/^{2~P)\ — 1 and p > 1, where a G (0,1/2) 

is to be determined. It is straightforward to see that for Sk = 2~P'^ + 

<2gV«(2-P.), 


To estimate the variation of Y on the intervals [sk, Sfc+i], define 


bk-.= inf r{t), := sup £{t), m*. := (a^-h 6 fc)/ 2 . 

te[sk,sk+i] te[sfc,si,+i] 


Note that r,i, and thus also /, are almost surely /3-Holder continuous for any {3 G (0,1/2) by, say, 
the Hoffman-Weilandt inequality (24). Let be the Holder constant for r,£. Assuming we choose 
(3 > a, there exists some random k£^r,i 3 so that for k < ki^r,i 3 , the path property f{sk) > g{sk) = {skY 
holds along with the estimate 1/2" — 2 c/ 3 |sfc — Sk+iY~°‘ > > 0, for some Then for such 

k 4 


2 c/ 3 |sfc|^ > /(Sfc) >bk-ak> f{sk) - sup |r(sfc) - r(f)| - sup \£(t) - £{sk)\ 

> g{sk) - 2c/3|sfc - Sk+if 

> gi2~PY -‘^c/slsk - Sk+iY (30) 

~ 20 "S/c-l-ll 2C/3|Sfc Sfc-|_i| 

^ Ca,l3\^k S/j-i-il 


Now write Ti '.= Sk and for j > 1 define stopping times 


aj := inf{f G [r,-, Sk+i] :W{t) = ruk}, Tj+i := inf{f G [aj, s^+i] : \W{t) - mk\ > {bk - afc)/4}. 

Let P^’'’, E^’'" denote the probability and expectation conditional on £, r and note B{t) is independent 
of £,r. Notice for all J > ^Tj+i]y = 0 (by Definition |2.1|) and we can estimate 

E'-’- [\W(c,) - M/(tj)|] < E'"- - “^0 + E'’’- 

< < St+l)[/(S|i+l) + /(St)]. 

Hence by Doob’s L^-submartingale inequality conditional on 3Ft , for all J> 1 , 


E^’^ [|y(a,) - y(r,)|] < E^’*- [|IH(a,) - W(r,)|] + E^’^ 


sup \B{t) - B{Tj)\l{ri<s^,+^) 


< [/(sfc+i) + f{sk) + 2{sk+i - ■ P^’''(rj < Sk+i). 
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Using now the Doob L^/"-subniartingale tail bound, we expect over the period [sk, only hnitely 
many Brownian oscillations greater than {bk — ak)/4: > Ca,/ 3 \sk — Sfc+i|“/4 (here, we used estimate 

S): 


Li>2 


i>2 j>2 


li.r 


sup \B{t) - B{aj-i)\ > 


Cq:,/ 3 I ^fc+1 I 


i>2 

E,>2(U 


(u - 


Ija 


C‘a,li\Sk ■Sfc+ll 

CT,_i 


< 


< 


(Cq,/^) ^ I'Sfc "Sfc-i-il 
41/" 

= : Q,r,a,/3 < 00, a.S. for k < ki^r,p, 




where the third inequality uses the assumption 1/(20;) > 1 with tj — aj-i < 1. Since Y is monotone 
on [u,(jj], j > 1 (by Dehnition B. our work so far yields, for k < ke^r,i 3 , 


E^’^’Vi 




F < [\Y(Oj) - F(t, 


i>i 


^ {Ce,r,a,l3 + 1)[/('SA:+i) + /(Sfc) + 2|Sfc — 

< {C£,r,a,l3 + l)[2C/3(sf_^^ + sf) + 2^’^^/^(2 

< {Ce,r,a,^ + 1)[2 ■ 4^c^(2-^’'=)^ + 2^+^/\2-P^y] 

< C {n-Vk\P 

by the /3-Holder continuity of / and the facts /(O) = 0 and c/ (2o) > c> (3. Summing this estimate 
over the relevant intervals [sfc/,Sfe/+i] contained in [2“^'=, 2“^'=+^] gives 


fc':K,,Sfc, + l]C[2-Pfc.2-Pfe+2] 


2-Pfc(2-Pfe)/^ 
dsfc'.Sfc'+i]-'' — '^e,r,a,p gl/a(^2-Pk'j 




Hence, as long as we take a,f3,c so that 1/2 < c < a{(3 + l),o < /3 < 1/2 (e.g., take o = 6/15,/3 = 
7/15), summing over p = Pk > ^ yields a hnite number, as required. 

To summarize, by repeatedly applying Proposition 2.5| and the conclusions of this proof, the 
pathwise constructed system (W(f))i>o of (14) is the unique strong solution to the non-degenerate 
interlaced system 


of reflecting Brownian motions with a wall started at the origin. Further, for 
all fc > 1, the kth. level is distributed as Pq’'^, each pair (W^’*'^, W^) distributed as and 

conditional on W^, the dynamics are independent of W^,..., for k >2. Convergence 

of symplectic Plancherel growth X{t,N) to this system is the content of Proposition 2.3, which 
completes the proof of Theorem |1.1[ □ 


3 Representation theory of symplectic groups 

3.1 Notations and ancillary results 

For n > 1, we write := [(n -|- 1)/2J. Let denote the set of partitions A = (Ai > ... > A^^ > 0) 
of nonnegative integers of length at most (the length of a partition A is the number of nonzero 
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terms). For A G Jn,/^ £ Jn+i, write A -< /i if A interlaces /i in the sense 


hi > Ai > /r2 > A2 > .. . > jam ^ > hr„+i- 

i arises naturally and frequently in representation 


For A G Jn, the transformation A* := Aj + r, 
theoretic formulas. For the construction of Section 4.1 
sequences u = (m^, ■ 


dehne the collection J, 


,u 


’') of partitions G Jj, 1 < i < n, of length n. Let J, 


n^seq 
n,paths C ' 


un,seg 

the subset of Gelfand-Tsetlin patterns, i.e., hnite interlaced sequences m = -< ■ ■ ■ -< m”). 

Joo.seq, ^00,paths denote inhuite versions of each of these sets. Recall our notational convention: 


of all hnite 
denote 
Let 


Important Notational Convention. Quantities with indices that overhow (e.g. A^^+i for A G J„) 
are set to 0 and those with indices that underhow (e.g. Aq for A G Jn) are set to 00 . 

For each r = consistently hx a basis {e_i, Ci,..., e_r, e^} for C^®' and dehne an inner product 
by [x, y] := X]fc=i ^kV-k — X-kVk for x,y E The even symplectic group Sp{2r, C) is the set 

of linear transformations T preserving this form: [Tx,Ty] = [x,y] for all x,y E C^®®. The odd 
symplectic group Sp{2r — 1, C) is dehned as the closed subgroup of transformations stabilizing e^, 
and we realize Sp{2r — 2, C) as the further subgroup stabilizing the pair e_r, e^. This furnishes 
a hxed choice of embeddings Sp{2r — 2,C) C Sp{2r — 1,C) C Sp{2r,C). Proposition 1.1 of |45] 
shows Sp{2r — 1, C) admits the (Jordan) semidirect product decomposition iL 2 r-i(C) x Sp{2r — 2, C) 
into unipotent and semisimple parts, where H 2 r-i{C) is the complex Heisenberg group. Dehne also 
the even compact symplectic group Sp{2r) := U{2r) fl Sp{2r,C) and dehne Sp{2r — 1) to be the 
subgroup of Sp{2r) stabilizing e^. 


Remark 3.1. An intrinsic, though diherent, dehnition for odd symplectic group is the set of 
transformations preserving a skew-symmetric bilinear form of maximal possible rank] see section 
9 of Proctor [H] for a discussion of candidate dehnitions. We have adopted the dehnition of 
Gelfand-Zelevinsky [2S] to exploit the above inclusions and to account for intermediate levels of the 
Gelfand-Tsetlin scheme (see [13] and Section 3.2). 


If n = 2r, the hnite-dimensional irreducible representations Vx of Sp{n, C) are parametrized 
by A G J 2 r- In the less common odd case n = 2r — 1, the same partitions A G J 2 r-i = ^ 2 r now 
parametrize certain indecomposable Sp{2r — l,C)-modules Z1(A) that exhaust the nonisomorphic 
representations of Sp{2r — 1,C) (see |13]). For either parity of n and for A G Jin, denote the 
corresponding character by Xni nnd write dim 2 rA := dim 5 p( 2 r.,c)hA, dim 2 r-iA := dim 5 p( 2 r-i,c)'^('^)- 


Section 3.2 shows these characters are expressible in terms the Ghebyshev polynomials of the second 
and third kind Jk,i/ 2 {.x), Jk-i/ 2 {,x), k > 0, which are dehned by the initial conditions Jo, 1/2 = 
1, Ji, 1 / 2 ( 2 :) = 2x and Jo,- 1/2 = 1, Ji,- 1 / 2 ( 2 :) = 2x — 1, while both satisfy the same three term 
recurrence ^ ^ 

xpk{x) = -pk+i{x) + -pk-i{x), k>l. 


(31) 


They also satisfy, for z on the unit circle. 


J, 


k,l /2 


z + z-^ 


^ —(/c+l) 


(32) 


Our choice of notation comes from the relationships 

T _ (2fc -F 2)!! ( 1 / 2 , 1 / 2 ) j {2k)\\ (- 1 / 2 , 1 / 2 ) 

2-(2A: + 1)!! ' {2k - 1)\\ ' 
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where are the Jacobi polynomials of parameter {a,/3) dehned to be orthogonal with respect 

to the weight W(a^p){x) := (1 — x)"(l + a:)^l[_i 4 ](a:). We also have = (^’^“) for fc > 1 and 

= 1 for fc = 0, so that Jk-i/ 2 {^) = 1 for A; > 0, a fact we will need to exploit many times. Lastly, 
dehne an inner product with respect to a normalized weight 


{f,9)a := 


20+1/2 


TT 


f{x)g{x)w(a,i/ 2 ){x)dx. 


(33) 


Using the relations in (32), one can compute directly Jk,a)a = which explains the choice 
of normalization in (33). We will frequently make use of the orthogonal decomposition 


T{X) = 

k=0 


(34) 


for T G C^[—1,1] (cf. Lemma 2 of [H]). See Szego [H] for more details on this discussion. 

We now collect some results that are fundamental for computing the correlation kernel in Section 
1^ and for uncovering the intertwining relationship in Section 

Lemma 3.2. The Chebyshev polynomials satisfy the following identities 
Yfr=Q Jr -1/2 = Js, 1/2, foT all S > 0. 

EZl‘^Jr, 1/2 = for all s > 1. 

2. For any T G C^[—1,1], 

j:Zs+l{Jr,-l/2,T).,/2 = (J.,1/2,T(1) -T)_V2 
YlT=s(Jr,l/2,T)l/2 = {Js-1/2, T)-1/2 

Proof. For x G [—1,1], let x = {z-\- z~^)/2 for some 2 ; on the unit circle. The hrst part then follows 
from using the identities in ( [3^ and exploiting the resulting hnite geometric sums. For the second 
part, the orthogonality relations and the hrst part together give 


{JsP/2, l)-l/2 — ('^s,l/2, Jo-l/2)-l/2 — ( Jr-1/2, Jo-1/2 ) — (Jq-1/2, J 0 -I/ 2 )-1/2 — 1- 

- 1/2 


I r=0 


Orthogonal decomposition (34) gives 

OO OO 

(Js,l/2, J'(l))-l/2 = J'(l) = ^(Jr,-l/2, J')-l/2Jr,-l/2(l) = ^ (Jr,-1/2, J")-1/2, 


r=0 


r=0 


where we have used the fact J^- 1 / 2 ( 1 ) = 1. Then subtract {Js,1/2,T)-1/2 = Yfr=o{Jr-1/2,T)-1/2 
from both sides to prove the hrst identity of the second part. For the other identity, note 

y^(Jr,l/2,T)i/2 = 2 { ^ r ^ 1-\ J") 1/2 = (1, T) _i/2 

r=0 ' ' 1/2 


where the term disappears necessarily from convergence of orthogonal decompositions (34) and 
the last equality uses the normalization of our inner products (33). Subtracting the hnite sum 
YfT=ii{Jr,i/2,T)\/2 = (1 - Js,- i/2,T)-1/2 from the right side completes the proof. 

□ 
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Define 


1/2, if n even 
— 1/2, if n odd 

For smooth E G C'“[—1,1], denote the mth Taylor Remainder of E abont 1 by 

E{x), m < 0 

E(x) - EETo' - 1 )‘. ™ > 1 ' 

and let := {x - l)’'"-'Rf 

Proposition 3.3. The functions n > 1, / G Z, satisfy the composition rule 




i{s,t) : = 


2-l(^<t), if n even 
l(s<i), if n odd 




t>o 


Proof. Fix n odd so r„_i = — 1 and = —1/2. Then nsing part (3) of Lemma 3.2 with 

Tm{x) := [x — l)~'^Rff{x) and the fact Tm(l) := lim^,.,.! Tm(x) = E^'^\l)/m\, we get 

OO 


t>0 


t=s+l 


= ( T 


4,1/2, 


Tr^-l Tr^_f(l) 

X — 1 


1/2 


= J. 


s,l/2, 




l-Tn-l 


(a: - 1) 




1/2 


where the second eqnality uses our inner product normalization (33). The case n even instead 


involves = r„_i, = 1/2, and part (4) of Lemma 3.2 


= 5 ^{a,i/ 2 ,(i - lY’-'RtJin 

r=s 

= (4_v2, (o; - l)'-"--'i?f_.„_,)-i/2 = K:\-iis). 


r>0 


□ 


We collect one last result that will help us express our probability measures on partition se¬ 
quences in determinantal form. 

Lemma 3.4. Indication of interlacement between two partitions A G Jn,/^ ^ Jn+i takes the deter¬ 
minantal form 

1(A^^) = det[0„(A,,/Ij)]["L\ X 

where Aj := Aj -|- — i. 


if n even 
if n odd ’ 
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Proof. First suppose \ jjL If u is even, then = r^+i — 1 while if n is odd, = r„+i. Under 
the transformation X, fi ^ X, Jl, the interlacing condition becomes 

j Pi > Ai > /I 2 > A 2 > ... > Pr-„+i > Ar„+i ^ -1, if n even 

\pi > Ai > p 2 > A 2 > • • • > Pr„+i > Xr„+i >0, if n odd 

Since = 2 ■ l(s<i) for n even, the strict inequalities p, > Aj imply cfniXijJlj) = 2 for 

1 < j < * < 'f'n+i) and similarly if n odd. Hence, the r„+i x r„+i matrix [0n(Ai, Pj)][^^A\ is triangular 
hlled with 2 ’s if n even and with I’s if n odd, which proves the statement. 

If A 7 ^ p, then let k* be the largest index of A such that the interlacing condition fails, i.e., one 

of Pa:*+i > Afc* or Afc. > p*,* holds. If A^* > p^*, then Aj > p^* > pj for 1 > f > fc* > j > r^+i, 

which implies that the block of [0n(Ai, with bottom left corner (A;*, k*) is hlled with O’s. If 

k* = Tn+i, then the k* column is a zero vector and ensures the determinant is 0. Otherwise, by 
maximality of k*, the k* and A;* + 1 columns are the same (the matrix is triangular), and so again 
the determinant is 0. The case Pfc*+i > A^* is handled the same way. □ 


3.2 Branching rules and character formulas 


Branching rules determine the decomposition of a representation when restricted to the action of 
a subgroup. But restricting an irreducible representation V\, X G J 2 r, of Sp{2r, C) to the action of 
Sp{2r — 1 , C) does not readily yield a decomposition, at least not into irreducible representations. 
Fortunately, the restricted V\ does decompose into a semi-direct sum of the indecomposable Sp{2r — 
1, C)-modules {£(p)}^gj 2 ^_j. More precisely, Shtepin |13] establishes the following (multiplicity-free) 
branching rules: for any irreducible representation Vx, X G J 2 r, of Sp{2r, C) and any indecomposable 
representation A2(A), A G ^ 2 r-i, of Sp{2r — 1, C), we have 


|sp(2r—1,C) —lA^(h)) A2(A) |5p(2r—2,C) ® 


(35) 


where ffl denotes semidirect sum. 

Fix A G Jn and recall the notation Aj := A, -|- — i. Proposition 5.1 of |13] shows that the 

indecomposable odd character X 2 r-i depends only on its semisimple part Sp{2r — 2, C). Using this 
fact with the branching rules (35) gives 


X2r-1 - X2r-ll5p(2r-2,C) - 


r-2- 


AtSj27 —2 


Note also the determinantal identities 


(36) 


det [z. 


i-1+1 .^-(*-1+1) 


][j=i = det[(zj + ^)* ■ Y[{zk - % ^) 


k=l 


det [z 


i-1+1/2 (i-l+l/2) 

“T x;.- 


][j=i = det[(zj + ■ W{zl'^ + 


1/2 -l/ 2 x 


k=l 


(37) 


Combining (32), (36), (37) with the Weyl character formula for the even case (see Chapter 24, [25]). 
we can compute for either parity of n and for M G Sp{n) 


det 


xtiM) = Xni^l, ■ ■ ■ , Zr) = 


J- 






*5=1 


det[(^j 


(38) 
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where the spectrum {zi, ... ,Zr„, z~^} of M has Zi on the unit circle and z^^ = z~^ = 1 if n 
odd. Finally, only for the sake of explicitness, the dimension formulas can be obtained either from 
the branching rules ([^ or from evaluating \imsp(n) 5 g^eXni 9 )^ where e is the identity: 


dim 2 rA 


n 


2 2 
mf — mi 

^ j 


n 

l<i<r 



dim 2 r-iA 


n 

l< 2 <j<r 


I? 


n 


m: 


/2 


m 


/2 ’ 


(39) 


where /j := Aj + (r„ — i) + 1 , := (r„ — z) + 1 , and l[ := U — 1 / 2 , m' := rrii — 1/2 for 1 < z < r„. 


3.3 Plancherel measures for Sp{oo) 


For a general topological group G, a function / : G —)• C is a class function if f{hgh~^) = f{g) 
for all g,h E G, and is positive definite if the matrix [f{gJ^gi)]2j=i is positive dehnite for any 
G G. Consider the inductive limit Sp{oo) := U'fLiSp{2r). A general characterization 
problem of Okounkov-Olshanski (Theorem 5.2, |39]; see also [19]) shows that the extreme points of 
the convex set of positive dehnite class functions y : Sp{oo) —)■ C normalized to have y(e) = 1 at 
the identity e are parametrized by the collection of triples (a, j3, 5) with a, [3 nonnegative decreasing 
summable sequences satisfying + ^ and 7 := 5 — + /dj) > 0. By slight abuse, we 

express these triples with either parameter <5 or 7 . For uj := {a, fi, 7 ) and M G Sp{n) with spectrum 
{zi, zf^,..., Zr„, zf^} (again, the Zi are on the unit circle and Zr^ = zf^ = 1 if zz odd), we have 

r„ 



where 




z + z 


-1 


= 2) n (i + f(^-i))(i + f(^--1)) 

ii(i-f(z-i))(i-f(z-i-i)) 


(40) 


(the odd case follows by observing the inclusion Sp{2r — 1) C Sp{2r)). Since the irreducible 
characters {x2r}x&s 2 ^ complete orthonormal basis for class functions on Sp{2r), we may 

write the restriction x‘^\n ■= x‘^\sp{n) as a convex combination 


X |2r = 


AGj 2 r 


xi 

dim 2 r A ’ 


(41) 


■) 2 r for characters. Evaluating both 


PfirW = 1 and positive dehniteness of x^ ensures that > 0. The P^ are 


where convergence holds with respect to the inner product 
sides at e yields X]AeJ 2 . 

thus probability measures on J 2 r- In the odd case, we cannot a priori rely on any classical theory for 
restrictions y^| 2 r-i to Sp{2r — 1), so instead we use (41) and the consequence (36) of the branching 
rules to get 

V- PnW 

rX2rl2r-l — 2^ 


X^\2r-l — (X^br) 


2r-l 


dim' 


2r 


AiSj2T--l 


AGj 2 r 

Pa").(A) dim 2 r -i /z 
dim 2 r A 


AGJ 2 ? 


dim 2 r A ^ 

11 ^ X 


X2r-1 


X2r-1 
dim2r-i fi 




P^r-M 


xi-1 

dim2r-i /i’ 


(42) 


The implicitly dehned Pfir-i are of course positive, and we will deduce unit mass on J 2 r-i directly 
from their explicit expression, which we derive next. For each u, we call the series Pfi, n > 1, the 
Plancherel measures for the infinite-dimensional symplectic group Sp{oo) and we refer to the Pfir-n 
r > 1 , as intermediate. 
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Theorem 3.5. Choose a triple u = (a,/3,7) to determine E‘^ as in (40). The associated series of 
Plancherel measures Pif, n>l, for Sp{oo) is explicitly given by 


P“(A) = 2('?) . det (x - 




dim„ A, A G 


(43) 


where ( 2 ) = 0 and we continue the notation Xj := Aj + — j for A G 

Proof. First assume n is even. Note that G C^[—1,1] and write Xi := {zi + z~'^)/2 G [—1,1], 
where zi,... ,Zr„ provide the eigenvalues for M G Spin). By orthogonal decomposition (34) and an 
inhnite extension of the Cauchy-Binet Formula (cf. Lemma 1 of [H]), we have 


det [{Xj - !)• , = det 




fc >0 


*J=1 


= 5^ det|{y (x - . det[Jfe,„.(i,)]5.,. 

A:i>...>A:r>0 

Making the substitution ki = fli, which takes a nonstrict partition fii > ... > frr„ > 0 to the strict 
one in the last summation above, and rearranging yields 

x“i„(A/)=n E“(x,) =i; det {X - 


k=l 


E det (x - ■ 


/iEJrz 


det [x 


r—tim 

j 


2 O det (x - ■ xMei..... zX 


/iEJrz 


where we have used det[(xj — 1)''"“®]]’"^! = det[xY~^YCj=i- Applying the inner product {■,Xn)n for 
characters to both sides and recalling PifiX) := {x‘^\n, Xn)n ' proves the result for n even. 

To compute for n — 1 odd, we use the formula just derived for Pifi^X) to get (note 

ttn-i = -1/2 and = r„_i) 


- E ^ = E det (- - 1)^"-^^")./. 


dim„_i /i dim„ A 


E det [(7,1/2. (^-1 )"-E")i/2 


J i,j=l 

1 Tn 

Pn-l 

5^(A,i/2,(a:-l)’'"--‘£“)i/2 

k=Jlj 

= det [(7,- 1 / 2 , (x - 1 )’'--B")_i/ 2 ];":\ , 


(/4-<A) 


det 


(/dj<Ai) 


AgJtj 


= det 


*J=i 


(44) 


*J=i 

rn-l 


where the third equality follows from our Lemma 3.4, the fourth from the generalized Cauchy-Binet, 
and the last equality follows from part (4) of Lemma 3.2[ 

For either parity of n, recall dim„/i = limsp(n) 3 g^eXn(,g) and write yj{g) := {zj{g) P Zj{g)~^)/2 
with the Zj{g) giving the eigenvalues of the limiting g G Sp{n). Using (38), we then conclude unit 
mass by computing 


E = E 2^’* ■ det [(7.«». (^ - . dim„/ii 


flG^n 


/i.EJn 
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lim 2, 


Sp{n)3g^e 


= lim 

Sp{n)Bg^e 


= lim 




det (a; - 1)’'" ■ det[J^,,„„( 2 /j(c/))][;;.^i 

det[fr'i9)]u=i 


det [EZo(Jk,a^^ - ir--^E-)^„Jk,M9))ZU 

det[y^-\g)]l-=i 

det [(yjig) - 


(45) 


Sp(n)3g^e det [(^/j (^f) - 1)^" 


lim Y[E^{yj{g)) = 1, 

5p(n)9g->-e -*--^ 

a=l 


where Cauchy-Binet and orthogonal decomposition (34) were used in the third and fourth equalities. 
This completes the proof. □ 


4 Determinantal correlation structure 


4.1 Central probability measures on partition paths 


The branching rules ( |35| ) imply that interlaced sequences u G Jn,paths ending in «”■ = A exactly 
parametrize the basis elements of either Vx or /1(A). We thus have the elementary, yet crucial, 
relation 

^ Wu = dim„ A. (46) 

'^^3n,paths 

u^=X 

where the weight Wu '■= l(ni-<a2_....^„n) = indicates interlacement. For u G Jn,seg, 

consider the cylinder sets Cu '■= {t G Joo.seglt^ = ,... = «"■}. We now realize the Plancherel 


measures on of Theorem 3.5 as embedded in a single probability measure on Joo,seq by 
the prescription 

^ •“ dim„n- ’ 

The weight Wu ensures P‘^ is supported on Joo,paths C Joo,seg- The consistency relations 

T’n-i(/i) 


dim„_i fi 


2-^ rH 


ACJn 


dim„ A 


(47) 


established as in (44) guarantee P'^ is well-dehned. Lastly given tt, n G Jn.seg, we obviously have the 
centrality condition P‘^{Cu)wy = P‘^{Cy)wu. 

The mapping 


^ oo,seg 


9 M !-)■ V{u) := = (■u™,m) | m > 1, I < k < r^}- 


pushes P^ forward to a probability measure := P^ oV~^ on which determines a simple 

point process in Z>o x Z>o, i.e., a random element of We call this process simply the 

Plancherel point process, though this title does not depend on the choice of coordinates. In these 
coordinates, the interlacing condition becomes 


X 


X 


n+1 

fc+1 

n+1 

fc+1 


<x1< if n even 

<x1< x'l'^^, if n odd 


1 < A; < r„. 
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and symplectic Plancherel growth can be visualized as 


3^3 

T-S 

X 2 Xl 

ks 

xl 

—>■ X 

xj 

xf 

^2 


xf 

3^2 

— 

-V 1 


xf 

Xl 

-t \X2 


xj 


1 

xf 


x\ 

-)■ 

1 

xf 



X 3 

X 2 


Xn 


X 2 x^ 

xf 


X, 


Xi 

xl 


Xn 


Xn 


Xn 


5 

X 2 

xf 


X, 


( 48 ) 




Define the kth correlation function ior zi,..., Zk G Zi>o x Z>o by 

^(^1, ...,Zk)-.= ^^{{E e |E D {^1,..., z,}}) = D {^1,..., z,}). 

The coordinatization of Figure is similarly determined by the mapping 

Soo,seq 9 M H- V{u) := {{y^ , m) = {uff + m-2k + l,m) | m > 1,, 1 < A; < r^}, 

which as before determines a simple point process in Z>o x Z>o via the probability measure 
:= o 'P~^ on 2 ^>oxz>o_ j|;g correlation functions pf are defined the same way. 

4.2 Computation of correlation kernel 

Theorem |1.21 of the introduction is a consequence of the following general statement. 


Theorem 4.1. Choose a triple u = (a,/3,7) as in Section 3.3 so that there exists a loop in C 


around the interval [—1,1] containing no zeros of the E‘^ of (40) (this is true if (di < 1). Then 
the correlation functions {pf}k>i of the general Plancherel point process are determinantal: for 

Zi,... ,Zk G Zi>o X Z>o, 

^{zi,...,Zk) :=P(T‘^ D {zi,...,Zk}) =det[K^{zi,Zj)]lj^^, 
where the (nonsymmetric) kernel is given explicitly by 

2 au+ll 2 1 rl r (rn\ 

K-{{s,n),{t,m)) = -— / f w I ^udx 

TT im J_iJ E‘^[u) [l — uy”^[x — u) 


oon+l/S rl 

+ l(rx>m)- / Js,aAx)Jt,ar.{x){l “ 

J-1 


(1 + xY^'^dx. 

(49) 

for {s,n), {t,m) G Z>o x Z>o, where the complex integral is a positively oriented (i.e., counterclock¬ 
wise) simple loop around [—1,1] containing no zeros of E‘^. 

For any n > 1, a point configuration X„ := {x™ GZ|l<m<?7,, 1 < k < r^} in Z>o x [n] 
determi nes the cylinder set := with u G Jn,seg satisfying = x^, 1 < m < n. Using 

our push-forward measure o V~^ takes the determinantal form 


3.4 


Lemma 

r (X„) = P^{Cu) = = 2^") n ■ det 






m=l 


dimn*^ 

/ n—1 

const n det | • det 


*J=1 


(60) 


km=l 
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where := {Js,a„^ (x —1)'’" ^Ri-r„)a„ was defined in Proposition 3.3 Also, define convolution 

over Z>o by (/ * g){x,y) := Y.z>o for bivariate functions f,g and by (/ * g){x) : = 

^z>o /(^> ^)9i^) if 9 is univariate. 

Proposition 4.2. For any {s,n), {t,m) G Z>o x Z>o cmd A; G Z, define the functions 


■■= F- 




2ni J E{w){w — l)^m-fc+i 


dw 




27ri 


J. 


S.OLn 5 


- 1 )’ 
X — u 


n I m 


du, for n < m, 


where the contours are positively oriented (i.e., counterclockwise) simple loops around [—1,1] con¬ 
taining no zeroes of E. Then the simple point process determined by of (50) has determinantal 
correlation functions pf with kernel 




(51) 


k=l 


Proof. The result is a restatement of Proposition A.2 so it suffices to identify the functions used 
in Appendix 1^ with the functions in the statement above. By orthogonal decomposition (34), we 
have for 1 < /c, / < 


5 ; _,(s) = f-i (x - «,.(») 

s>0 s>0 


E{w){w — Ifixi-i+i 


dw 


27ri J {w — 1)*^ 


dw = 6ki. 


Note also that is a polynomial in t of the same degree as 1, t), where 

Qk = 2A: — 1. These items conhrm that {<h™ -ddYkZi i® unique basis of the linear span of 
1 , A)}fc=i that is biorthogonal to the 

We also need to show = 0 ^ * ... * cfrn-i for n < m. First assume m = n + 1. If n is odd, 

rn = rn+i, an = -1/2, On+i = 1/2, and 4>n{s,t) = l(s<t), so that 


27ii 


■Js-1/2, 


Jtp/2{u) 


X — u 


) du — {Js-l/2i Jt,1/2)-1/2 
/ -1/2 


"A,—1/2) ^ ^ A', —1/2 j 0n('5) t)i 

-1/2 


r=0 


where the hrst equality computes the residue at u = x, the second follows from (1) of Lemma 3.2 


and the last from orthogonality. Similarly, if n is even, findfi) = 2 ■ l(s<t), = f'n+i — 1) and the 

a.’s switch, so that 



i-l 


2 / t/s,l/2) ^ ^ A-,1/2 J t)i 

\ r -=0 / 1/2 

where now the hrst equality computes residues at both u = x and u = 1, and uses the fact that 
Jt- 1 / 2 Z) = 1- The full statement for general n < m then follows by induction. □ 
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From our current (51), a few more calculations are required to arrive at our main expression 


(49) for in Theorem 1.2 


Lemma 4.3. For any n, m > 1, 


{u - - E{u)) F{x- ly-E 


{x — u)E{u){u — 1)' 


J. 


(x-1) 


Fn Fm 


• m ' n 


S^Oin 5 


E 


Ji 


t^OLn 


Jt,am {u^du 


(52) 


Note how the terms in (52) simplify if > r^, since in this case = E. 


Proof. First note the identity 


E 

k=l 


u — 1 
X — 1 


U — 1 


X — U 


U — 1 
X — 1 


Then if r^<rn, 

Fm 


k=l 


J. 


SiO(.m 


2711 J E{u) 


J. 


S^OLn 1 


{x — 1)^" E{x) 
{u — 1 )'’™ X — u 


u — 1 
X — 1 


(53) 


du. (54) 


and taking the residue at n = a; of the second summand in the inner product yields (52). If instead 
rm> rn, write 

Fm Fn Fm 


k=l 


k=l 


k=r„+l 


The hrst summand is treated as in (54): 


1 


k=l 


2ni J E{u) 


Js,a^ / . {x- ly- E{x) 


{u — !)'■"* X — u 


1 - 


u — 1 
X — 1 


du. 


For the second summand, we use the identity ( |53| ) repeatedly to get 

Fm / 1 \ A* 

z;r ~ J-(r^>r„+l) 

k=rn+l ' ^ fc=r „+2 r=l 

'u — 1 


= (E(x) - E(l)) 


(x — 1)^-'’" 

(u - 1)^"+1 


^(rm>rn+l) 2 -^ 2 ^ 


r\ {x - ly-'^^-^ 


1 - 


X — U 

Fm—Fn — 1 


- 1 


{rm>rn + l) 


E 

r=l 


X — 1^ 

E ^^ yi ) (m-i)^"+'-+i 


1 - 


X — U 


and then arrive at 




fc=r„+l 


277i 


\ ds,ani 


{x — 1)’ 


Fn Fm 


' m ' n 


U — 1 
X — 1 


du 


Tm-rn-r ' 


1 

271 i 


\ J.‘ 


E{u){x — u) 

{u - iy--^^{R^^_,yu) - E{u) + E{x)) 


du 


E{u){x — u) 

Taking the residue at u = a; of the hrst term and combining with the hrst summand completes the 
proof. □ 


































Now plug the expression (52) into (51). If the expression (52) simplihes, quickly leading 


to the main expression (49) for in Theorem 1.2 (note we have multiplied by the conjugating 
factor “(— 1 )'’"-^™”^ which vanishes in the determinant). But if we get (49) along with the 

additional term 


1 / / . . . ^ -r (W) \ , / r r 


m ' n 


The residue of the hrst term at u = a; exactly cancels the second, and since 


k=rm-rn 


k\ 


there is no residue at u = 1 (recall by assumption E = E'^ does not have a pole at 1 either). This 
completes the proof of Theorem EH 


5 Multilevel Markov process 

We state and prove the results of this section without coordinatizing Joo,seq 


5.1 Single level dynamics 

The point processes arising from our general construction in Section are parametrized by 
the triples oj = (a,/?, 7 ) of Section 3.3 The next proposition dictates a natural way to transition 
between the Plancherel measures on a single level Define by replacing E'^ with 'ip in the 
expression for of Theorem 3.5 For any fixed p) G C^[—1,1], define a x J„-matrix with 
entries 

^ dim„ A 


T*(i,, X) := det 






dim„/i' 


Proposition 5.1. Let be nonzero at 1. 


A) = P^'^^iX), and similarly 

2 - = 

3. If ipm G 1,1] converge uniformly to if, then T^™(/i, A) —>■ X) as m ^ 00 . 

4 . Tfp has nonnegative entries if ipfx) = po + Pix with po > pi > 0 or if ijj^{x) = 7 > 0 . 

5. The (stochastic) semigroup operating on the Banach space of absolutely summable 

functions /^(Jn) on J„ is Feller: \\t¥ — 0. 


For example, set uq = (0, 0, 0), and for /c > 1, let Uk be given by a = (g, g,..., g, 0,...) (i.e., 
k q’s), (3 = (0, 0,...,), and 7 = 0 . Since E‘^'‘ ■ E^^ = P‘^'=+i for fc > 0, this choice along with the 
proposition yields a discrete-time Markov chain X*, t > 0 with state space and distribution Ppf^ 
at time k. 
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Proof. Both items of the hrst point follow readily from the generalized Cauchy-Binet and orthogonal 
decomposition (34); for example, 






= 2 ( 2 ) dim„ A ^ det ({x - 1)'’" > 1 ) 


■ det 


= 2 ( 2 ) dim„ A ■ det 
= 2 ( 2 ) dim„ A ■ det 






J i,j=l 


’ 00 > 


^ k=0 




'^A„a>2, (a; - 1)’'" >1 




*j=i 


The second point is similar to (45): use dim„ A = limg^eXn(5') (38) to compute 




lim 


AEJti 


Sp{n)Bg^e 

dim„ fi 


det ( 


det [ 4 . (%(c/) )]["=! 


AEJn 


ij=i det[(2;j(5() + 


-i Wn-* 


Ji J = 1 


_ hm 5 p(„) 9 g_,.e det Efc=o(d/i*,an' 0 ) '^fc,a„)Q:„'/fc(|/j(5'))]j j = i 
dim„p det[(zj(c/) + 

dim^p det[(z,(^) + ^ 

where are the eigenvalues of the limiting g G Sp{n), with yj{g) := {zj{g) + Zj{g)~^)/2, and 
where Cauchy-Binet and orthogonal decomposition ( [3^ were used in the second and third equalities. 
The third point of the proposition is a straightforward application of the dominated convergence 
theorem. 


Now for the fourth point, hrst take 'i/j{x) = po + pix. Using the three term recurrence (31), 
orthogonality relations, and xJo,a„ = + Jo,a„{^ — 2 Q;n)/ 2 ], we compute 


Pi 

{Jjli,ar,-, J\^^ajp)ar, = P0^Jli,Xj + 




(1 - 2 an) 




(55) 


If Xi > Pi+l for some 1 < ^ < r„, then A^ > pi + l ior 1 < k < i < I < Vn, which implies T^{p, A) = 0 
(the resulting matrix admits a 2 x 2 block form with an off-diagonal block of O’s and a diagonal 
block with a zero vector). The same conclusion of course holds if pi > Xi + 1, so assume |Aj —/ij| < 1 
for all 1 < i < r„. If \Jli — Jli+i\ > 1 for some 1 < i < r„, then |/Ij — Aj+i| < 1 implies |Aj — Jli+i\ > 1 
and similarly |Aj —/Ii+i| < 1 implies \Jli — Xi+i\ > 1. In either case, T^{p, A) breaks into a product of 
determinants, one of which is the i x i upper left corner minor of [(T^i,a„, T 5 ;, Q,^'0)a„]ij=i- Iterating 
this argument reduces consideration to the case where \ — < 1 for all 1 < i < r„, which means 

Pi are all equal to some p G Z>o. Now the two blocks corresponding to {i : Aj = p±l} are triangular 


with nonnegative entries (55) and straddle a tridiagonal block corresponding to {i : Aj = p = pi]. 
Thus, writing q* := max{z : Xi = p = pi}, := min{i : Xi = p = pi}, q := q* — q* + 1, we have 
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further reduced consideration to the determinant of the q x q tridiagonal matrix 


A(A) := 


Po pi/2 0 0 

Pi/2 Po pi/2 0 

0 pi/2 Po pi/2 

0 . 

0 . 


0 

0 

0 

: 0 
Po Pi/2 

Pl/2 Po+P 2 ,g- 


(56) 


where p 2 ,q* ■= ^ p. Writing r± := {po ± ^/pl — p\)/2 and hn{x) := pox^ — ^x^ we can 

compute, using the tridiagonal determinant recurrence relations and initial conditions, 


det Ag(A) = 




pI 


[hq{r+) - hq{r_) + p 2 ,r (r+ - r”)] . 


(57) 


Notice that h„(^) = 0 and h is increasing on [|^,cx)). Since r+ > r_ > we have hq{r+) — 
hq{r_) > 0 and so the expression (56) is positive, as reqnired. The case now follows 

from the previous parts by noting 


T'^~i = lim = lini ... 2 "i+ 7 ( 3 ;-i)/m_ 


( 68 ) 


For the hnal point, note that stochasticity of implies 


luy - /<i||,.(J») = sup E - ry(A, fi)) = sup(2 - 2Ty(A, A)) 


AgJ- 


/iEJn 


AGJn 


and (57) implies 


inf T^+l(x-l)/m(^ ^ hn{r+) hn(r_) 

A 6 J„ " ^ ^ 


I -Pi 


n—1 \ 2 / 77—2 

k \ I \ ^ ^(n—2)—k \ 7^0 

4 


=po\Yl ^+^- 


,k=0 


E 

. fc =0 


r_^rl 


'4“l, 


( 59 ) 


where we have taken po = l — ^/m>'^/m = pi in dehnitions that depend on these qnantities. By 
(58), Tn^(A, A) > A)]™, so letting 7 —f 0 in (59) completes the proof. □ 


5.2 Intertwined multilevel dynamics: discrete steps 

Dehne cotransition probabilities from J„+i down to level by 

dim„ A 


rr‘(p.A)^ .. " 


cf. the consistency relations (47). Identity (46) implies dim„+i/i is the nnmber of interlaced 


sequences of length n + 1 ending in p, so summing T”+^(/r, A) over A G Jn evaluates to 1 . These 
stochastic operators ensure that the interlacing condition is preserved when we link single level 
dynamics to a mnltilevel evolntion, bnt the following intertwining relationship is key. 
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Proposition 5.2. Fix ifj G C'^([— 1,1]) with '0(1) = 1- For n>l, the stochastic operators and 
satisfy the intertwining relations 

A 71+1_ rTin-\-l rjilp _ rpll) rTin-\-l 


Proof. First assume n is even, so that = 1/2, On+i = —1/2, Tn = r„+i —1, and (j)n{s,t) = 2-l(s<t). 
For pi G Jn+i and A G Jn, we compute directly 




•2GJn 

dim„ A 
dim„+i /i 


Z £ J fi 


dim„ 
dim„+i /i 


■i(«,)rt(z.A) 


2 ^det|.#>„(j,,/^i)|J"+\det 


Z^^n 


+ , 1 / 2 ) +^-, 1 / 2 ^’ 


1/2 


J i,j=l 


where we have used Lemma 3.4 Expand the first determinant along the r„+i-column, where the 
convention ^r„+i = — 1 applies. Omitting constants, the Zth resulting summand, 1 < / < Tn+i, is 


^ ^ det /li)]i<i^i<rn+i det 


2GJtj 


^<j<Fn 


+i.l/2) 


1/2 


J i,j=l 


= det 


= det 


= det 


= det 


= det 


= det 


^ ^ 4 >nik) /ij) (jk,l/2y 


k>0 


1/2 


l<i^l<r„+i 

i<j<rn 


‘^Jk- 1 / 2 , +j,l/ 2 V^ 
fc =0 / 1/2 


^_1 >+,, 1 / 2 ^ 


l/2_ 






2 ^(+i,-l/2'0)(l) +i,-l/2'0) +j,l/2^ 


- 1/2 




2 ( +. _l/2'0, ^ Jk-1/2 

\ k=Xj-\-l ' _]^^2 

oo 

1/2^07 ^ ^ 07i("^_;7 1/2 

/c =0 


i<j<rn 


- 1 / 2 . 


l<i<?’n 


where we have used Cauchy-Binet in the first equality, the first part of Lemma |3.2 for the third, 
multilinearity/skew-symmetry of determinants and the assumption (+.^_i/ 20 )(l) = 1 for the fourth. 


and the second part of Lemma 3.2 for the hfth. Continuing our calculations, we sum over I to get 


dim„ A 
dim„+i p, 


Tn+l 


2 -r-n+i ^(_l)r„+i+Z 


1=1 


+^,-1/21/) 'y ^ 0n(Aj, k)Jk—1/2 


k=0 


- 1 / 2 . 


l<i<»'n 
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dim„ A 
dim„+i /i 

dim„ A 
dim^+i n 


det 


Jfii,— 1 / 2 '^1 ^ ^ 0n(Aj, k^Jk,—1/2 


k=0 


-l/2_ 


r„+i 


*J=1 


2 ^-+1 det (J^^.,-i/ 2 > 


-1/2 


zeJ, 


^n + l 


■ det 


n+1 


^J = l 


'^'i) 


n + 1 

J ij=l 


= (r„*+iT„"+‘)(f,,A), 


where for the hrst equality we note that the entries with j = r„+i are all 1, which follows from 


Ar„+i = —1, {Jji.-i/ 2 'ip){l) = 1, and orthogonal decomposition (34). The much more straightforward 

□ 


n odd case involves the other components of Lemmas 3.2 and 


3.4 and is left to the reader. 


For M, t G 


> n,paths 


dehne 




{u,t) : = 






(60) 


In words, it is the probability of the transition —)■ —)■ by a jump and cotransition 

conditional on —)■ occurring by such steps; note this quantity only depends on 

Let ip be either {l—pi)+pix, 1/2 > pi > 0, or := Just as the stochastic operators 

account for transitions between the probability measures on the stochastic (by Proposition 


5.1) transition operator 


n—1 




k=l 


supplies an evolution of probability measures on Sn,paths of the form 


n—1 


pLm,m pi(n ■ IT pPh ”*'’«e 


( 61 ) 


k=l 


in the usual sense {Ptlaths^t)^'^) = P^palhs ^> 1,^2 e <^^[-1,1], ^ 0. To prove this, 

compute, for t G Sn,paths, 


iPtathsAtm= Y 

paths 

n—1 


n—1 




k=2 

n—1 




k=l 




E !?■(«“) 

k=3 

/ k=2 

\ n—1 n—1 

Y P) n = ppp{p) n t^), (62) 


n—1 




k=l 


k=l 


where the last equality follows from the hrst part of Proposition 5.1 
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5.3 Second construction of symplectic Plancherel growth 

Proposition 5.3. Let and let ip G C^[—1,1] satisfy '0(1) = 1. Then 


pV’ plQn _ pi’ Ai’^y 

n^paths^ n.vaths^n ’ 


where the paths x ^n^paths-i^citrix Qn is the infinitesimal generator of (uncoordinatized) symplectic 
Plancherel growth considered up to level n. 

Proof. Let 53„ be the Banach space defined as the completion of the snbspace of l^{^n,paths) con¬ 
sisting of measnres of the form (61) corresponding to fnnctions nonzero at 1. The stochastic 
operators form a semigronp on by (62). The Felle r pr operty for follows 

from the same property of (the fifth part of Proposition 5.1): the form (61) implies 


E 

paths 


n—1 


pi’ii’2 ( \ 
^n^paths V ^) 




5^ E 


k+l(yk+l 






— paths ^—1 


\Pp*Hn'‘) - Ppiu")], 




SO we have 


\\Af - = r*-' - '’4° 0 


■'0 

n,paths 


^'yQn — 


Hence (see, e.g., Chapter 19, 0 ). there exists an operator Qn on 53^ snch that P, 

F’t,paths^^ 1 SO suffices to show Qn = = Qn- Since the calculation for the second order 

approximation of only differs from m for transitions involving jumps into the wall, we only 
present this case, the others being similar. 

Assume the system is in a state so that a wall jump is possible at an odd level A; > 1 , i.e., a state 
t e In,paths that satisfies = 0, > 1, so > 2, and if fc > 3, > 1 (cf. 

(48)). Consider a transition to u G Jn,paths that agrees with t G In,paths except that ut. = ut. =1. 
We compute 




(63) 


to second order. Noting fi-y{x) = 1 -|- 7 (x — 1) + 0 ( 7 ^) for small 7 > 0, we get, similarly to (56), 


T^(^^n^)=det 


'^t",-l/2) '^ap-1/2'0^ 


= det 


(1-7) 




2 p,t^+i 


2 t|+i ,0 


0 

0 


= 7-2-(^+i/2) + 0(72). 


- 1/2 

0 


pfc 


*j=i 


(1-7) 


2 V+i,t§ 


2 t^tf+i 
(1-7) 


0 

0 


2 tLtJ+i 


( 1 - 7 ) 

0 


0 

0 

0 

0 

0 


sy . 2-C+V2) 


0(7^ 


34 













where [3 = ±1/2 is the second parameter of Jacobi polynomials (Section |3.1[ ), indicating whether we 
are dealing with the initial conditions of the hrst (-) or third (±) kinds of Chebyshev polynomials, 
which correspond to the orthogonal and symplectic case, respectively. Similar calculations in the 


denominator of (63) yield an order of 1 ± 0 ( 7 ). Since the dimensions occurring in the numerator 
and denominator can all be shown to cancel, we arrive at 

7 ■ 2-(/^+V2) + 0 (j^) 




1 ± 0(7) 

= [7 ■ ± 0(7^)] (1 ± 0(7)) = 7 ■ ± 0(7^). 


where the second equality uses the geometric series. This shows that the orthogonal case /3 = —1/2 
involves rate 1 wall jumps and the symplectic case (3 = 1/2 involves rate 1/2 wall jumps, as 
required. □ 


6 Asymptotic analysis of symplectic Plancherel growth 

6.1 Bulk limit: incomplete beta kernel 

In this section, we work in the coordinates 9 A h->■ A. Focusing on the kernel K'^ at two limiting 
values, assume 7 > 0 and (s, n), (t, m) G Z>o x Z>o depend on N in such a way that 7 ~ Nt > 0, 
s, f ~ Nv > 0, and r„, ~ Nrj. Assume the differences s — f, Vn — rm are of constant order. Recall 

the double integral term from our correlation kernel K'^: 


20 , 1 + 1/2 


TT 


27r/ 






- Jsa {x)Jta {u) , 


{x — !)'■" (1 — x)""(l ± 


X — U 


dudx, 


where the positively oriented w-contour encloses the unit circle. Deforming this w-contour to be a 
loop centered at 1 and passing through cos^, for some 9 to be determined, will produce a residue 


that is considered below. With an aim toward using the identities (32) make the change of variables 
u = (n±x“^)/2 so that the x-contour is outside the unit circle but connecting to Similarly, 
with the change of variables x = [z + z~^)/2, the weight becomes 


(1 - l)“(l + xY'^dx m„{dz) := \ _ i /o 


(64) 


Siz 


The additional factor of “1/2” appears because the mapping is two to one. Taking into account the 
residue at x = x from the hrst deformation, we are left with the kernel 


ri] t, m) 

2 Q!„+ 1/2 


,7P+^-1)/2 


2 tiH 7 / 2 '^+"" 


z + z 


-1 


Ji 


t,(Xr} 


v + v ^ — 1)^" (1 — X '^)dvma{dz) 


^ v-\-v 


±1 


20n + l/2 


— 1)'’’" 2 : ± z~^ — (x ± X“^) 

( 66 ) 


(n>m)' 


TT 


J. 


|z|=l 


Z + Z 


-1 


Ji 




z + z 


-1 


z + z 


-1 


- 1 


' n • m 


ma{dz) 


2 o„+1/2 


71 


J. 


s,a„ 


Z + Z 


-1 


Ji 




± X 


-1 


z + z 


-1 


• n ' m 


(66) 


niaidz). 


where the last contour connects clockwise along the unit circle. 
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Lemma 6.1. The following limit holds: 


(■n — 


lim (66) = l(n>m)- 

N^OO '-' 


2 ni 


kl=i 

— Oim 


271 i 


^{t+rm,+oim.)-{s+rn+oi„)-l _ 2:'j'^~'^dz 

^{t+rm+am)-is+rn+an)-l _ 2 :^'^~'^dz 


where the contour connecting is negatively oriented and crosses (—oo, 0 ). 


Proof. Both terms in (66) can be handled by the same argument, so assume n > m and concentrate 


on the hrst of these. We will use the identities (32) but by shrinking or enlarging the unit circle, 


we may ignore terms involving This leaves us with 


20^71 + 1/2 


-J. 


71 


z + z 


-1 


Ji 


t.OLri 


z + z 


-1 




z + z 

2 

Tn-Tm 


-1 


ma{dz) 


= 




dz ^ 

Siz TT 
I’n— Pm. 


S — t\ ( Z+Z 


= < 


(^Zt-s ^ 

^^t-s+1/2 _ ^s-t- 1 / 2 ^ _ xj 

2 ;{t+rm)-{s+r„)-l ^x _ ^'^‘^'^ri-‘2rm 
_^(t+rm)-{s+rn)-‘S/2 ^x _ 


dz ^ 

A.iz TT 

Pn Pm 


if (y^i O^m 2 
d ^2 if Q, = -a =1^ 

21 / 2 + 2 - 1/2 Siz n 2 

if 0)2 O'm 2 

W!±£Zh!k£i if a =-o =-- 
2 — 2 : ^ 422 TT n m 2 


(67) 


dz 


2^n-rm27r2 
2Pn—‘^Pm 1 — 2 dz 


if CX^i Clm 2 




(l-z) 


2 rr? — 2 r 


n m 


•c _ _ 1 

^172 2 ’"rr-rm 27 ri ~ “ 2 

if otfi otjYi 2 


ti2 


y(i+7’m)-(s+r„)-l/2 


(1-z) 


2'-n-rm27ri 
2rn-2rm dz 


l-z 2^ri-rm2m 


if 0)7, 0)77, 


where the second equality follows hrst from making the substitution z —)■ z~^ to the integrals 
determined by the second summands in each case, then from using the identity + z~^ — 2 = 
z~^{z— 1)^. Note that in the hrst of these steps, the transformation changes the orientation of the z 
contour over the unit circle but an additional sign change occurs after reorienting. Letting N ^ 00 
completes the proof. □ 


For the double integral term (65), we again use the identities (32) to compute 

Js,a^ ( Jt,a^ 


( -dz 
822 
—dz 
Siz 


= < 


[z^'v^^zv — iz + z — z^v * 

^ ^ ^ 'Z ZV 


2 I z—z ^ 
tij v — v~^ 


s„,—t z 


= 


dz 

Aiz 

dz 

Aiz 

( -dz 
Siz 
—dz 
Siz 

dz 

Aiz 

dz 

Aiz 


{z'^v^)z^ — (z — Z ^^^ + Z^V 

(z'^v^)^/zv + (z~^v^)^ + Z~^V~^^= + 

^ ^ V \ J y/zV y/v 


{z''v^)^/zv + {z ^ 

-2(;2-V)^ + 2 z-^v-^^ 

-2(z-V)^ - 2 z-^v-^^ 

2(z-V)^ + 2z-"u-‘^ 

2(z-V)^ - 

'' ' y/Z yf ZV V 


— 1/2 

^l/ 2 ^^-l /2 

^l/ 2 +^-l /2 

^l/ 2 ^^-l /2 


l/2_^^-l/2 

2l/2+^-l/2 

^l/ 2 +^-l /2 

^l/2 + ^-l/2 


if c^ji 2 

if O^Yi 2 

if 2 

if CXt). O^rn. 


if 2 

if 2 

if 2 

if CX-n. OtrY), 


( 68 ) 
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where the second equality follows from the substitution ^ —)■ z~^ applied to the integrals determined 
by the hrst and fourth summands in each case. Now we make the additional transformation v —)■ v~^ 
applied to the integrals of the second terms in each case, which completes the n-contour into a 
positively oriented loop through the points as indicated by the dotted line in Figure Most 
importantly, we have a double integral over two loops with integrand 


2«n+l/2 / ^ 

2 nH 


-1 


Ji 


t,(y.n 


2 Jz + z^ — {v + v'^j 


1 2 ; 


1 — V 


-2 


2 : V—V~^ 
Jv z — z 


if O-n CTm 2 


2 m 2 'Kiz z + z 


-1 


(v + v 


-1 


-dvdz ■ < 


•r _ _ 1 

^„l/ 2 +.y-l /2 II <^n — ~Otm — 2 


n 

yl„l/2+.„-l/2 if ttn — «m — — 2 
,1/2+ -1/2 ^ ^ _i 

^ yjz V — V ^ 2 


The part of the integrand depending on the hydrodynamic scaling parameter N becomes 

^-N{S{zo)-S(z)) 


(,N{S{v)-S{zo)) ’ 

where Zq is to be determined and where, as in ([^, 


S{z) = := r 


z + z 


-1 


— v log z + r] log 


z + z 


-1 


-1 


(69) 


(70) 


To identify the correct steepest descent/ascent deformations, we hrst rely on Proposition 5.1 of [H] 
to determine 6 and zq\ 

Proposition 6.2. (Proposition 5.1 o.f IW) For t/t] >1), define 


q±{r,r]) = q± 




®i + i± 

V 





3/2 


but set to 0 if 1/2 < T/rj. Let Zmax,Zmin be the largest, smallest real roots of S\z), and let go be 
the complex root of S'{z) m El \ D . Also define 


{ jyfrozen | ^ . j/ > 7 ^ . g_|_(^x, 7 /) } 

jyhqmd _ | ^ ]^3 fj ■ q_{T,rj) < U < 1] ■ g+(r, T])] 

pempty | ^ ]^3 ; < 7 ^ . g_ (^x, 7/) } 


Then the quantity 


Zo = Zo(t, Z/,7/) 


fzmin, (T,iy,p) 

< go, (r, 12,7/) e 

(r,7/,7/) ePf™ 


is critical for S{z) and satisfies Zmax > 1 in the first case; |go| > 1 in the second; and Zmin < —1 in 
the last. 


Our work below will conhrm that the dehnitions of the regions j^frozen^ j^empty gjygj^ 

here in this proposition are consistent with the introduction’s. 
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Figure 6 : The yellow region signifies the region lH((S'(n) — S'(zo))) > 0 and the white elsewhere. 
The dotted blue line indicates the completion of the n-contour into a loop by the transformation 

V v~^. 


For (r, z/, rj) G choose 6 so that are points on the unit circle in the region of ^ where 

the condition 9 d((S'(; 2 o) ~>S'(^))) > 0 holds. Since |^o| > 1, we enlarge the ^-unit disk to be a steepest 
descent loop staying in the region —fH((S'(zo) ~ ^'(z))) < 0 but passing through the critical points 
zo and Similarly, make the steepest ascent n-loop pass through the points zq,^ and stay in the 
region !H((S'(n) — S'( 2 ;o))) > 0. Figure [^exemplifies these deformations. The double integral (65) is 
now seen to vanish as —>■ oo, but leaves behind residues. By inspection, the residue ai v = z~^ 
involves the term with |n| < 1 for Lebesgue almost all v and hence vanishes as iV —)■ oo. The 


expression (69) tells us the residues at v = z are 


^ r' .-.-1 

4 V 2 


1 if (Xfi 2 

if Q( = — (y. = — 

^ 11 LXfi 2 

1 if (X^i 2 

V 2—1 2 


i r.—Ti±id_i 

2TTi Lie \ 2 


' n • m 


1 if Oiji 2 

^ if = -ar,, = I 

1 if Otji CKm 2 

—-— if O' = —o = — — 

V 2—1 2 


• dz 


( 71 ) 


We have implicitly used the fact that —)■ 1 as n —)■ 2 ;. Hence, we have proven 


Lemma 6.3. //(r, z/, 77 ) G then 


n '-*-m 


lim (65) 

V-s-oo - 


+ 


2711 


y(t+rm+am)-is+rn.+an)-l 


J ZQ 

Oifi — 


27ii 


yit+rm+oim)-is+rn+a„)-l 
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Noting that the conjugating factor 1)“"““™-” vanishes in the determinant, Lemmas 

|6.1| and |6.3| complete the proof of Theorem 1^ once we conhrm the regions are consistently dehned 
with the introduction. 




Figure 7: The plot exemplihes steepest descent deformations for the frozen (left) and empty (right) 
regions. The yellow region again signihes the region yi{{S{v) — S{zo))) > 0 and the white elsewhere. 


For the case (r, i/,ri) G do not deform the initial u contour and only make the substi¬ 

tutions X = {z + z~^)/2, u = {v + v~^)/2, as before. Steepest descent deformations as in the plot 
of Figure do not produce residues ai v = z~^ since zo{t, v,rj) = Zmin < — 1 and also no residues 
at n = z by considering the indicated regions of the plot. Hence, the double contour vanishes as 
—)■ oo without leaving behind any terms and we are left with the triangular matrix governed by 
the hrst term in Lemma 16.II 


^n>m)' 


■1)^ 


2 Tii 


z 


(t+rm+Olm)-is+rn+an)-l 


(1 - z)^-^dz. 


Pl=l 


Since the diagonal entries of this triangular matrix are 1, the determinant of converges to 1. 
This conhrms the dehnition of 

Finally, for (r, i/,ri) G ^ again refrain from deforming the initial u contour, but make the 

substitutions x = {z + z~^)/2, u = {y + n“^)/2, as before. We may perform the same calculation 


leading to (68), but now the transformation n —>■ n ^ to the second summand at the end of (68) 


leads to a loop inside the unit circle (currently the z contour). This second summand disappears 
without leaving anything behind since we may shrink its n-loop arbitrarily to avoid residues at 
V = z~^ while making z-deformations as in Figure But these deformations in the hrst summand 
of (68) leave behind a full loop of residues at v = z similar to 0 and Lemma [6.3[ 




27ri 


yit+rm+am)-is+rn.+an)-l 


(1 


^dz, 


Pl=l 


Since we did not perform any deformations to the initial u contour, we are asymptotically left with 
the sum of this last residue and the hrst term in Lemma 16.11 


-1 


(ri<m )' 




2 TTi 


^{t+rm+cim)-(s+r„+an)-ln _ 


P|=l 


dz 
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Since the diagonal entries of this triangular matrix are 0, the determinant of converges to 0. 
This conhrms the dehnition of 


6.2 Edge limits 

6.2.1 Finite distance from the wall: Jacobi kernel 


Assume again that 7 ~ iV ■ r > 0 and r„, ~ fV ■ 7 > 0, but that s, t are hxed and hnite. Assume 

only the difference n — m is of constant order. Let A{z) := rz + ? 7 log(l — z) and note 1 — rj/r is a 
zero of A'{z). Write our kernel as 


2 “"+i /2 1 r^re-N{A{i-v/T)-A{x)) ^ lx - 1 )’’- (1 - 

n ^ J_^j e^Ain)-A{l-rj/r)) Js,aAx)Jt,aM dudx 


X — u 




Now deform the w-contour, as in Figure to be a steepest ascent loop remaining in the region 



Figure 8: The plot exemplihes steepest descent deformations when targeting a hnite distance from 
the wall. The yellow region signihes where 91 (A(m) — A(1 — p/r)) > 0 and the white elsewhere. 


7l{A{u) — A{1 — rj/r)) > 0 and passing through the critical point 1 — 7 /r. If 1 — r^/r < —1, the 
unit interval already lies in the region TZ{A{1 — rj/r) — A(x)) > 0, so the double integral term 
tends to zero without picking up residues. As for the frozen region above, the kernel converges to 
a triangular matrix whose diagonal entries are I’s. But ifl — 7 /r>—1, then we acquire a residue 
at M = X given by 


2an+l/2^_]^y 




TT 


- XY^ + xfl'^dx 


'-1 


Noting that (—1)'^” is a conjugating factor completes the proof of Theorem 1.4 
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6.2.2 Corner of phase transition: Pearcey kernel 

Assume that 7 ~ N/2, and s ~ > 0, t ~ > 0, and lastly that Vn — N ^ > 0, 

Tm - N ^ y/N T]2> 0 

Lemma 6.4. Let x = — 1 G [—1,1]. Then for large N, 


2“Vx' 


Proof. Write x = cos 9 and note the power series 

r 1 


cos ^{z — 1) = TT 


^/2 


— , - =dx = TT - y/2y/z + H-, 

0 v^Vl-^/2 12 


If a = -1/2, then by _i/2(a;) = so that 

,cos[(s + l/2)7r — + o(l)] 


(-l)V,,_v2(:r) = (-1) 

= (-l) 


cos[7r/2 — ^Py/x' 12 + o(l)] 

cos[(s + l/2)7r + 0(1)] cos[i/i\/^] + sin[(s + l/2)7r + o(l)] sin[z/i\/2^'] 


cos[7r/2 + 0(1)] cos[A^“^Ay/rY2] + sin[7r/2 + o(l)] sm.[N-^Py/x'/2] 
sin[viy/^] 


[N-P^^/Pj 2 ]' 


sin 


and we may similarly use Js,i/ 2 {.x) = if a = 1/2. Using limj^_5.o = 1 completes the 

proof. □ 

Recall for any u = {a, f3, 7), the complementary point process of is also determinantal 

with /cth correlation function 

^’^(^1, ...,Zk) :=F {{X^y D {zi,..., Zk}^ = det[A'X(^i, u)]*^j=i, 

where the kernel K'^ is readily computable (using orthogonality for the Chebyshev polynomials) as 

2“"+V2 1 r E^{x) ^ ^ Ax-iy-{l-xY-{l + xf/^ ^ , 

TT 2711 J_i J E^[u) [U — ly^ X — U 

- l(„>m) (1 - . 

(72) 

Deforming the u-contour in the double integral term of (72) (with uj = (0,0,7)) in Figure 
we endeavor to hnd the contribution at —1. Making the substitutions x' = N^^‘^{x + 1) and 
u' = + 1), we have, for large iV, 


(1 -a;)“"(l + _ ^^.^{2 - x' 


du'dx' 


X — u 


O^n .» /rpl 

^-3/4 ■ du'dx'. 


Putting this together with Lemma 6.4, we get 




X — u 


u'{x' — u') 
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Figure 9: Steepest descent deformations when targeting the corner of the phase transition. The 
yellow region signihes the region y{{A[u) — 74(—1)) > 0 and the white elsewhere. 


The remainder of the integrand satishes 

g 7 a; g-Af(/l(-l)-^(a;))+»;i\/7V(log(l-a:)-log2) l%-r)\x'(2 

X — 'f'n ^ •' _ _ . 


— 'f'nO'f'm — in 


^ (^ll — lym gAr(A(u)-A(-l))+772\/]V(log(l-M)-log2) g-(M')^/8-^2«72 ’ 

where A{z) := zj^ + log(l — z) and the last approximation follows from, respectively, second and 
hrst order Taylor expansions abont —1. Hence, pntting everything together, we have 


^1/42C 


■m — 


-2)' 


(- 1 ) 


t — S 


2«n + l/2 ^ pl 


71 2711 


'-1 






[x — lY" (1 — xY"{i + 


dudx 


-)■ 


^2 1 


200 POO , /.o , /v2 / / 

f e 8 2 ^ sin[i/i\/^] sin[z/2\/^] 


u — X 

dx'du' 


71 277% 

Tnrning to the other term of note the following two asymptotic relations: 


u'{u' — x') 


{-2Y^-^Yx-iy 


n ' m 


0 2^/iv) 


{vi-V2)VN 


-(Vl-V 2 )^' 
—>•6 2 


2"’"(-l)Vs,a„(a;)(-l)*Jt,a^(a:)(l - x)“"(l + xY^'^dx ~ N sin[viV2x'] sin[1/2 

yx' 

Then the single integral term satishes 


-)■ 


1 I f°° -(^ 1 -^ 2 )^' . r /lZ~n ■ \ 

\ -/ ^ ^ sm[z/iV2a;'J sm[z/2V2a:'J^^ 

' '0 vx' 


= 1 
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1 


= 1 


\/277 Jo 
1 


{vi>V2) 




-(Vl-V2)^' 

e 2 


exp 


COs[(z/i + 122) V^] — COs[(z/i — 1/2) \/^] 

(z/i + z/2)^i r {1Z1-1Z2Y 


dx' 


Vl - V 2 


exp 


hi -h2 
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where the last equality follows from standard Gaussian computations. Noting that the conjugating 
factors will disappear in the determinant completes the proof of Theorem 1.5[ 


A Algebraic argument for determinantal correlations 


For completeness, this section covers the linear algebraic details on which Proposition 4.2 relies; 
we follow a combination of Theorem 4.2 of [U] and Lemma 3.4 of m, which in turn rely on the 
Eynard-Mehta theorem in the manner of |l8j . 

Fix n > 1. For 1 < m < n and N G Z>o, consider point sets ;= of variables x™ 

ranging in {0,..., iV}, where we suppress dependence on N. Dehne ^ := and recall our 

work in Section |4.2| furnishes a measure on point conhgurations in up to level n by 


n—1 


const. n det ■ det [d-” ■ 


(73) 


m=l 


Note (73) uses our notational convention = — 1 for m even (recall the proof of Lemma 3.4), 


so add the virtual variables V := {vi = —as a notational device to form the larger space 
:= V U 1^. To separate these variables out, let hFm, 1 < m < n — 1, be the x matrix 

with entries xj^’*'^) := 0m(xj", xj^"*"^), and similarly, for 1 < m < n, dehne a V x matrix 

Em-i to have entries Em-i{vi, x'p) = a^y^)l{i=r™)l(m odd)- Dehne a x V matrix to 


have entries := “ 1)^'“ Proposition 

Now the measure (73) is proportional to a minor of the X x X matrix 


3.3 


L := 


0 Eq 0 F/2 

0 0 -ITi 0 

0 0 0 - 1^2 


Q 


En -1 

0 

0 


-W^-1 

0 


0 A 
C D 


The matrix L supplies a point process (called an L-ensemble) on the larger space X with distribution 
V e 2^ where Lx is the X x X submatrix and Ix '■= Ixxx- Ds correlation 

function is determinantal with kernel Kx = Ix ~ {Ix + L)~^ (see [ 22 ]). Write I^ + L 

where as indicated, A is the Vxblock, C is the x V block, and D the ^x^ block. One can see by 
inspection that the 2) x matrix D~^ has (i, j)-block W[ij) := (hFi ■ ■ ■ hF,_i)l(i<j) + IxU)xx(r)'^i^=j)■ 
'Notice 

/ 2 3 n '' 

AD ^ = ( i7ofF[i^i), Ek-lW[k, 2 ), Ek-lW[k,3), • • • r Ek-lWlk,n) 


k=l 


k=l 


k=l 


Dehne also an x matrix Mn := AD 


Proposition A.l. The point process of particle configurations up to level n governed by (73) is 
the same as the if) conditional L-ensemble defined by the prescription Pg)(D) := det[i^+L^] ’ ^ ^ 
which is determinantal with correlation kernel 


K = L 


2 ) 


{hj + L) bx2) - L 


9 


D 


-1 


+ D-^CM-^AD 


-1 


(74) 
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More explicitly, its x -block is given by 


^{ij) ■= 


M-^ I 


1 >^J) 


,k=l 


(75) 


Proof. For the the basic facts on conditional ensembles and in particular the hrst equality of (74), see 
. The second equality is the content of the Eynard-Mehta theorem, and is proved by computing 


the inverse of the 2x2 block expression of Jg + L = 


0 A 
C D 


in order to compute (Jg + L) ^Isgxsy- □ 


Now let iV —)■ cx) so that matrix multiplication induces convolution over Zi>o, which is defined by 
(/ * g) (x, y) := E^>0 y) for bivariate functions /, g and by (/ * g) (x) := J 2 z>o ^)g(^) 

if g is univariate. We need to compute (75). Let := W[i^n) and note := Ixo). Then for 
each 1 < j < n, 

' j \ j 

'^Ek-iWikj) ) {Vs,xl) = k odd){(pk-i*(p^’'’^^){vk,xl) = (0g,_i (76) 


.fc=i 


k=l 


where qk ■= 2k — 1. Take a V x matrix with entries denoted ^^{vs,xl) =: ixl)li s<rj), 

whose nonzero rows determine a basis for the row space {(0g^_i * a;)| x E of 

and which is biorthogonal to the in the sense 


E = is,, l<m<n. 


k >0 


Define an Vj x rj, 1 < j < n, matrix to have entries 

Bit ■= (0g«-i * l<s,t< rj. 

Add O’s to extend to an x matrix B^, which acts as a change of basis in the sense 

B^^^ = [f)o * 02 * ..., 0 , 0 ,...]' 

(note we identify n* EE i). Note further that for 1 < f < r^. 


x>0 


where the second equality uses (by Proposition 3.3) and the last uses orthogonal 

decomposition and -1/2(1) = 1. This implies in particular that if B^~^ is upper triangular for 
j > 2, then so is BP Hence, B^ is upper triangular and since B^ = B^ = Mn by dehnition, we 
have 


M, 


-1 


Ek-iW[k, 


J) 


,k=l 


B^ * 
0 * 


-1 


B^^^ = 


\B 


^■)-i *'■ 
0 *' 


B^ 0 
0 0 




We have thus computed the (i, j)-block (75) of the kernel to be 




M-^ [YEk-iW^k,j) 


,fc=i 




(See [TT] to hud a slightly more complicated expression for the kernel when the upper triangularity 
of Mn fails.) The entrywise statement reads as follows. 
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Proposition A.2. For any {s,i), {t,j) G Z>o x the point process determined by (73) has 
determinantal correlations with kernel 


K{{s,i), {t,j)) 




k=l 
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